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I. INTRODUCTION 

The realization of Bose Einstein condensates and 
quantum degenerate Fermi gases with cold atoms has 
been one of the highlights of experimental atomic 
physics during the last decade 1 -, and in view of re- 
cent progress in preparing cold molecules we expect 
a similarly spectacular development for molecular 

cn „ cm ]^ ca 2, 3, 4. 5, 6, 7, 8, 9 , 10, 11, 12, 13, 14, 15, 16, 17, 18. 19, 20, 21, 22, 23 

The outstanding features of the physics of cold atomic 
and molecular gases are the microscopic knowledge of 
the many-body Hamiltonians, as realized in the exper- 
iments, combined with the possibility to control and 
tune system parameters via external fields. Examples 
are the trapping of atoms and molecules with magnetic, 
electric and optical traps, allowing for the formation of 
quantum gases in ID, 2D and 3D geometries, and the 
tuning of contact inter-particle interactions by varying 
the scattering length via Feshbach resonance o 24 ' 25 . This 
control is the key for the experimental realization of 
fundamental quantum phases, as illustrated by the 
BEC-BCS crossover in atomic Fermi gase a 26 i 27 ' 28 i 29 ' 30 , 
the Kosterlitz-Thouless transition 3 - 1 - and the superfluid- 
Mott insulator quantum phase transition with cold 
bosonic atoms in an optical lattic o 32 ' 33 . A recent 
highlight has been the realization of a degenerate 
magnetic dipolar gas of 52 Cr atom a 34 ' 35 i 36 . Below we 
will be mainly interested in heteronuclear molecules 
prepared in their electronic and vibrational ground 
state. The new feature of polar molecules is their 
large electric dipole moments associated with rotational 
excitations. The new aspects in a condensed matter 
theory of cold polar molecules are the large dipole-dipole 
interactions between molecules. This points towards the 
possibility of manipulating these strong and long-range 
interactions with external DC and AC microwave fields. 
In particular, this raises interesting questions of cold 
ensembles of polar molecules as strongly correlated sys- 

f crrl :: 37, 38. 39. 40.41. 42.43.44.45. 46, 47, 48. 49. 50. 51. 52. 53. 54. 55. 56. 57, 58. 59 

Magnetic or electric dipole moments in external fields 
can give rise to anisotropic, long-range dipole-dipole in- 
teractions. Analogous to the case of cold atoms with con- 
tact interactions, a key element for the realization of in- 
teresting quantum phases and phase transitions with in- 
teracting dipolar gases is the capability of controlling and 
tuning system parameters using external fields. Much 
work has been recently devoted to the study of cold 
collisions in dipolar gase o 60 ' 61 '^ 2 -^ 3 -'^'^'^, which in this 
book is reviewed in the chapters contributed by Hut- 



son, Bohn, and Dalgarno. In the context of degener- 
ate molecular gases a significant body of recent work has 
focused on the regime of weak interactions, where the 
isotropic contact interaction potential competes with the 
anisotropic long range dipole-dipole interaction. For ex- 
ample, the existence of rotons in weakly-interacting dipo- 
lar gases has been predicted^^^^iIlSiSiZiiS, while 
exciting prospects have been envisioned for rotating sys- 
tem a 76 ' 77 ' 78 ' 79 ' 80 ' 81 and polar molecules in optical lat- 
ticea 82 1 83 J 84 1 85 J 86 1 87 1 88 1 89_ Below we will be mainly in- 
terested in the many-body dynamics of polar molecules 
in the strongly interacting limit. In particular, we will 
discuss a toolbox for engineering interesting many-body 
Hamiltonians based on the manipulation of the electric 
dipole moments with external DC and AC fields, and 
thus of the molecular interactions. This forms the basis 
for the realization of novel quantum phases in these sys- 
tems. Our emphasis will be on condensed matter aspects, 
while we refer to the contribution by Yelin, DeMillc and 
Cote in the present book for application in the context 
of quantum information processing. 

This review is organized as follows. In Sect. [IT] wc give 
a qualitative tour through some of the key ideas of en- 
gineering Hamiltonians and of the associated quantum 
phases. This is followed by two slightly more technical 
sections, Sect. IIII Al and Sect. IIIIB1 where we provide 
some details of the realization of a 2D setup where par- 
ticles interact via purely repulsive 1/r 3 potentials, and 
where we sketch how to design more complicated inter- 
actions by using a combination of AC and DC fields. 
Finally, Sect. IIVI deals with the applications of the engi- 
neering of interaction potentials in the context of the re- 
alization of strongly correlated phases and quantum sim- 
ulations. 



II. OVERVIEW: STRONGLY INTERACTING 
SYSTEMS OF COLD POLAR MOLECULES 



In this section we give a qualitative overview of many- 
body physics of cold polar molecules with emphasis on 
strongly interacting systems. In the following sections 
we will return to the various topics in a more in-depth 
discussion. 
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FIG. 1: (a) System setup: Polar molecules are trapped in the (x, t/)-plane by an optical lattice made of two counter-propagating 
laser beams with wavevectors ±1cl = ±fcLe z , (blue arrows). The dipoles are aligned in the z-direction by a DC electric field 
Edc = EocBz (red arrow). An AC microwave field is indicated (green arrow). Inset: Definition of polar (i?) and azimuthal (tp) 
angles for the relative orientation of the inter-molecular collision axis ri2 with respect to a space-fixed frame, with axis along 
z. (b) Qualitative sketch of effective 2D potentials Vjp(p) for polar molecules confined in a 2D (pancake) geometry. Here, 
p = ri2 sin #(cos tp, sin tp, 0) is the 2D coordinate in the plane z = and p = ri2sin$ (see Fig. [T] inset). Solid line: Repulsive 
dipolar potential V^p(p) = D/p 3 induced by a DC electric field. Dash-dotted line: "Step-like" potential induced by a single 
AC microwave field and a weak DC field. Dashed line: Attractive potential induced by the combination of several AC fields 
and a weak DC field. The potentials V^g (p) and the separation p are given in arbitrary units. 



A. Effective many body Hamiltonians 

Hamiltonians underlying condensed matter physics of 
N structureless bosonic or fermionic particles have the 
generic form 



i=l 



2m 



+ K 3 ff D ({r,}), (1) 



where pf/2m is the kinetic energy term, and Vtrap(ri) ls a 
confining potential for the particles. The term V^f ({r^}) 
represents an effective vV-body interaction, which can be 
expanded as a sum of two-body, three-body interactions 
etc., 



N 



N 



V*? ({*}) - Yl V™ (v i -r j )+ W 3B (r. t , r„r fc ) + . . . , 

i<j i<j<k 

(2) 

where in most cases only two-body interactions are con- 
sidered. Typically, V^ D ({r^}) must be interpreted as ef- 
fective interactions valid in a low energy theory, obtained 
by integrating out the high energy degrees of freedom of 
the system. In the following we will discuss the deriva- 
tion of such many body Hamiltonians for a cold ensem- 
ble of polar molecules in the electronic and vibrational 
ground state, and show how the effective interactions in 
these systems can be "designed" on an vV-body level via 
control of rotational excitations with external fields, in a 
form which is unique to polar molecules. 



Our starting point is the Hamiltonian for a gas of cold 
heteronuclear molecules prepared in their electronic and 
vibrational ground-state, 



N r 



2m 



TV 



£Vdd(ri-rj). 



(3) 



Here the first term in the single particle Hamiltonian cor- 
responds to the kinetic energies of the molecules, while 
Vtrap(ri) represents a possible trapping potential, as pro- 
vided, for example, by an optical lattice, or an electric or 
magnetic traps. The term describes the internal low 
energy excitations of the molecule, which for a molecule 
with a closed electronic shell : S(^ = 0) (e.g. of the type 
SrO, RbCs or LiCs), correspond to the rotational degree 
of freedom of the molecular axis. This term is well de- 
scribed by a rigid rotor = H^ t = BJf with B the 
rotational constant (in the few to tens of GHz regime) 
and J; the dimcnsionless angular momentum. The rota- 
tional eigenstates \ J,M) for a quantization axis z, and 
with eigenenergies BJ(J + 1) can coupled by a static 
(DC) or microwave (AC) field E via the electric dipole 
moment d^, which is typically of the order of a few Debye. 
For distances outside of the molecular core, the two-body 
interaction as given by the second line of Eq. is the 
dipolc-dipole interaction 



V dd (r) 



d l ■ dj - 3 (dj • e r ) (e r • dj) 



(4) 
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Here r = |r| = — r^] denotes the distance between 
two polar molecules, with e r the unit vector along the 
collision axis. In view of the large electric dipole mo- 
ments this term provides a comparatively strong, long- 
range anisotropic interaction between the molecules. 

The many body dynamics of cold polar molecules is 
thus governed by an interplay between dressing and ma- 
nipulating the rotational states with DC and AC fields, 
and strong dipolc-dipole interactions. In the absence of 
electric fields, the molecules prepared in a ground rota- 
tional state J = have no net dipole moment, and inter- 
act via a van-der-Waals attraction V v dw ~ ^Ce/r 6 , rem- 
iniscent of the interactions of cold Alkali atoms in their 
electronic ground-states. Electric fields admix excited 
rotational states and induce static or oscillating dipolcs, 
which will interact via strong dipole-dipole interactions 
Vdd with the characteristic 1/r 3 dependence of Eq. ([?]). 
Note that two parallel dipoles repel each other, while 
dipoles aligned along the collision axis will attract each 
other, possibly inducing instabilities in a many-body sys- 
tem. Thus to obtain stable many-particle phases is often 
possible only in reduced geometries, i.e. in combination 
with an external trapping potential Vtrap(ri). Finally, 
we emphasize that microwave excited rotational states of 
polar molecules arc long-lived, which makes these states 
available without the penalty of introducing significant 
decohcrence. This is in contrast to atomic systems, where 
spontaneous emission from laser excited electronic states 
is one of the main contributions to decohcrence. 

The connection between the full molecular A^-particle 
Hamiltonian (|3|) including rotational excitations and 
dressing fields, and the effective Hamiltonian |T]) can 
be made in a Born-Oppcnhcimcr approximation. If we 
diagonalize for frozen spatial positions {r^} of the N 



molecules the Hamiltonian^ Hbo = 

J2i<jV<id( r i — r j)> we obtain a set of energy eigenval- 
ues ({i"i}), which are interpreted as the effective N- 
particlc potential in the single channel many-body Hamil- 
tonian (fTJ). The dependence of V^f on the electric 
fields E provides the basis of the engineering of the many 
body interactions of two-body, three-body terms etc. in 
@. The validity of this adiabatic approximation and 
of the associated decoupling of the Born-Oppcnhcimcr 
channels will be discussed below. 

The above considerations set the stage for a discus- 
sion of engineering many-body Hamiltonians for polar 
molecules, and associated quantum phases. In the fol- 
lowing subsections wc will discuss specific examples of 
DC and AC field configurations for designing two-body 
and three-body interactions, which will be the content of 
Sects. Ill Bl and III CI below. Our discussion can also be 
adapted to optical lattices, and thus to a derivation of 
Hubbard models for polar molecules (Sect. Ill D[) . Fur- 
thermore, we can extend these derivations to molecules 
with internal degrees of freedom representing a spin 
fSect. lHE]) . Extended Hubbard models with couplings to 
phonons arc presented in Sect. Ill Fl for molecules trapped 



in self-assembled dipolar cystals. 



B. Self-assembled crystals. 

The conceptually simplest example, although remark- 
ably rich from a physics point of view, is a system of cold 
polar molecules in a DC electric field under conditions of 
strong transverse confinement. The setup is illustrated in 
Fig. [lja). A weak DC field along the z-direction induces 
a dipole moment d in the ground state of each molecule. 
These groundstate molecules interact via the effective 
dipolc-dipole interaction V^P(r) = D(r 2 — 3z 2 )/r 5 ac- 
cording to their induced dipolcs, with D = d 2 . For 
molecules confined to the x, y-plane perpendicular to the 
electric field this interaction is purely repulsive. For 
molecules displaced by z > the interaction be- 

comes attractive, indicating an instability in the many 
body system. This instability is suppressed by a suffi- 
ciently strong 2D confinement with potential a X / trap (z i ) 
along z, due to, for example, an optical potential induced 
by an off-resonant light field^. 

The 2D dynamics in this pancake configuration is de- 
scribed by the Hamiltonian 



rr2D _ 

-"off ~ 



2m 



(5) 



which is obtained by integrating out the fast z-motion. 
Equation is the sum of the 2D kinetic energy in the 
x,y-planc and a repulsive 2D dipolar interaction 



V£(p)=D/p* 



(6) 



with = {xj—Xi,yj—yi) a vector in the x, y-plane [solid 
line in Fig. [TJb)]. The distinguishing feature of the sys- 
tem described by the Hamiltonian ([5]) is that tuning the 
induced dipole moment d drives the system from a weakly 
interacting gas (a 2D superfluid in the case of bosons), to 
a crystalline phase in the limit of strong repulsive dipolc- 
dipole interactions. This transition and the crystalline 
phase have no analog in the familiar atomic bose gases 
with short range interactions modelled by a pseudopo- 
tential of a given scattering length. 

A crystalline phase corresponds to the limit of strong 
repulsion where particles undergo small oscillations 
around their equilibrium positions, which is a result 
of the balance between the repulsive long-range dipolc- 
dipole forces and an additional (weak) confining potential 
in the x, y-planc. The relevant parameter is 



_ -Epot 



E 



kin 



D/a 3 
h 2 /ma 2 



Dm 



(7) 



which is the ratio of the the interaction energy and the 
kinetic energy at the mean interparticle distance a. This 
parameter is tunable as a function of d from small 
to large. A crystal will form for 1, when in- 

teractions dominate. For a dipolar crystal, this is the 
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FIG. 2: (a) Tentative phase diagram in the T — rd plane: crystalline phase for interactions r«j > r QM and temperatures below 
the classical melting temperature T m (dashed line). The superfiuid phase appears below the upper bound T < Tvh 2 n/2rn 
(dotted line). The quantum melting transition is studied at fixed temperature T = 0.014D/a 3 with interactions = 5 — 30 
(dash-dotted line). The crossover to the unstable regime for small replusion and finite confinement u± is indicated (hatched 
region), (b) Strengths of the dominant three-body interactions Wijk appearing in the Hubbard model of Eq. @ for different 
lattice geometries: (bl) one-dimensional setup; (b2) two-dimensional square lattice; (b3) two-dimensional honeycomb lattice. 
The characteristic energy scale Wo = 72-D-Ro/a 6 is discussed following Eq. (|41|l . 



limit of large densities, where typically collisions become 
harmful. However the crystalline phase will protect a 
cold ensemble of polar molecules from (harmful) close- 
encounter collisions. This density dependence is in con- 
trast to Wigner crystals with 1/r- Coulomb interactions, 
as realized e.g. with laser cooled trapped ionsSi. In this 
case r c = (e 2 / a) / ti 2 / ma 2 ~ a and the crystal forms at 
low densities. In addition, the charge e is a fixed quanti- 
tiy, while d can be varied as a function of the DC field. 

Fig.[2ja) shows a tentative phase diagram for a dipolar 
gas of bosonic molecules in 2D as a function of and 
temperature T. In the limit of weak interactions < 1, 
the ground state is a superfiuid (SF) representing a finite 
(quasi-)condensate. In the opposite limit of strong inter- 
actions rd 3> 1 the polar molecules are in a crystalline 
phase for temperatures T < T m with T m « 0.09-D/a 3 , see 
RefjSl. The configuration with minimal energy is a trian- 
gular lattice with excitations given by acoustic phonons. 
In Rcfi^ we investigated the intermediate strongly inter- 
acting regime with > 1, using a recently developed 
Path-Integral Monte-Carlo technique (PIMC)^, and we 
determined the critical interaction strength r QM for the 
quantum melting transition from the crystal into the su- 
perfiuid. We found the latter to occur at r QM = 18 ± 4, a 
result which has been confirmed with a number of quan- 
tum Monte-Carlo techniques (author?)—'—. 

In Sect. IIVI below we return to a more detailed dis- 
cussion of these quantum phases, and in particular of 
the crystalline phase, and we show that the relevant pa- 
rameter regime where these phases occur is accessible 
with polar molecules. Besides the fundamental interest 
in dipolar quantum gases, the crystalline phase has in- 
teresting applications, e.g. in the context of quantum 



information^. We will return to self-assembled dipolar 
lattices below in a discussion of Hubbard models. 



C. Blue-shielding and three-body interactions 

By combining DC and AC fields to dress the manifold 
of rotational excitations we can design effective interac- 
tion potentials V 3D (ij — r 3 -) with (essentially) any given 
shape as a function of distance. For example, the ad- 
dition of a single linearly-polarized AC field to the con- 
figuration of Fig. [3a) leads to the realization of the 2D 
"step-like" potential of Fig. [TJb), where the character 
of the repulsive potential varies considerably in a small 
region of space. The derivation of this effective 2D inter- 
action is sketched in Fig. [3] and it is discussed in more 
detail in Sect. IIIIB2I below^^. The (weak) DC-field 
splits the first-excited rotational (J = l)-manifold of 
each molecule by an amount M, while a linearly po- 
larized AC-field with Rabi frequency O is blue-detuned 
from the (\g) — |e))-transition by HA, see Fig. E][a). Be- 
cause of hS and the choice of polarization, for distances 
p (cP/M) 1 / 3 the relevant single-particle states for the 
two-body interaction reduce to the states \g) and |e) 
of each molecule. Figure GD[b) shows that the dipole- 
dipolc interaction splits the excited state manifold of the 
two-body rotational spectrum, making the detuning A 
position- dependent. As a consequence, the combined en- 
ergies of the bare groundstate of the two-particle spec- 
trum and of a microwave photon become degenerate to 
the energy of a (symmetric) excited state at a character- 
istic resonant (Condon) point pc = (d 2 /hA) 1 ^ 3 , which 
is represented by an arrow in Fig. [3£b). At this Con- 
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FIG. 3: Design of the step-like potential of Fig. [TJb): (a) Rotational spectrum of a molecule in a weak DC field. The DC field 
splits the (J = I )-manifold by an amount hS. The linearly-polarized microwave transition with detuning A and Rabi frequency 
O is shown as an arrow, (b) BO-potentials for the internal states for £7 = (bare potentials), where \g;e}± = (\g; e) ± |e; g))/y/2. 
The resonant Condon point pc is indicated by an arrow, (c) AC-field-dressed BO-potentials. The dressed grounstate potential 
has the largest energy. 



don point, an avoided crossing occurs in a field- dressed 
picture, and the new (dressed) groundstatc potential in- 
herits the character of the bare ground and excited po- 
tentials for distances p 3> pc an d p "C pCi respectively. 
Consistently, Fig. [3fc) shows that the dressed ground- 
state potential (which has the largest energy) is almost 
flat for p ^> pc and it is strongly repulsive as 1/p 3 for 
p "C pc-, which corresponds to the realization of the step- 
like potential of Fig. QJb). We remark that, due to the 
choice of polarization, this strong repulsion is present 
only in the plane z — 0, while for z / the ground- 
state potential can turn into attractive. Thus, here the 
optical confinement along z of Fig. [lja) is necessary to 
ensure the stability of the collisional setup. 

The interactions in the presence of a single AC field 
are studied in quite detail in Refi££, where it is shown 
that in the absence of external confinement, this case is 
analogous to the (3D) optical blue-shielding developed in 
the context of ultracold collisions of neutral atom a 96 ' 97 ' 98 , 
however with the advantage of the long lifetime of the ex- 
cited rotational states of the molecules, as opposed to the 
electronic states of cold atoms. The strong inelastic losses 
observed in 3D collisions with cold atom a 96 ' 97 ' 98 can be 
avoided via a judicious choice of the field's polarization, 
eventually combined with a tight confinement to a 2D 
geometry (as e.g. for the case of Fig. [3] above). For exam- 
ple, in Refj^ 9 - it is shown that for interactions in the pres- 
ence of a DC field and of a circularly polarized AC field 
the attractive time-averaged interaction due to the rotat- 
ing ( AC-induccd) dipolc moments of the molecules allows 



for the cancelation of the total dipolc-dipole interaction. 
The residual interactions remaining after this cancelation 
are 'purely repulsive 3D interactions with a characteristic 
van-der-Waals behavior V c 3 ff D (r) - (d 4 / HA) / r 6 . This 3D 
repulsion provides for a shielding of the inner part of 
the interaction potential and thus it will strongly sup- 
press inelastic collisions in experiments. This will pos- 
sibly lead to the realization of quantum degenerate sys- 
tems of molecules and maybe of tightly-packed crystalline 
structures in three- dimensions^. 

A cancelation of the leading effective two-body in- 
teraction similar to the one sketched above in a dense 
cloud of molecules can lead to the realization of systems 
where the effective three-body interaction W 3U (ri, Tj, r/.) 
of Eq. (|2|) dominates over the two-body term V 3D (rj— rj) 
and determines the groundstate properties of the sys- 
tem. This is interesting, since model Hamiltonians with 
strong three-body and many-body interactions have at- 
tracted a lot of interest in the search for microscopic 
Hamiltonians exhibiting exotic ground state properties. 
Well known examples arc the fractional quantum Hall 
states described by the Pfaffian wave functions which ap- 
pear as ground states of a Hamiltonian with three-body 
interaction s] 100 ' 101 ' 102 . These topological phases admit 
anyonic excitations with non-abelian braiding statistic. 
Three-body interactions are also an essential ingredient 
for systems with a low energy degeneracy characterized 
by string nets ; 103 ' 104 which play an important role in 
models for non-abelian topological phases. This possi- 
bility of realizing a Hamiltonian where the two-body in- 
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teraction can be manipulated independently of the three- 
body term has been studied in Reb^ 9 -. There, it is shown 
that a stable system where particles interact via purely 
repulsive three-body interactions can be realized by com- 
bining the setup above with the tight optical confinement 
provided by an optical lattice. In fact, the latter serves 
the two-fold purpose of ensuring the collisional stability 
of the setup and of defining a characteristic length scale 
(the lattice spacing) where the exact cancelation of the 
two-body term occurs. The details of this derivation are 
given in Sect. [TVCl below, in connection with the deriva- 
tion of an extended Hubbard model with three-body in- 
teractions introduced in the following Sect. Ill Dl 



D. Hubbard lattice models 

Hubbard Hamiltonians arc model Hamiltonians de- 
scribing the low-energy physics of interacting fermionic 
and bosonic particles in a lattice 1 ^ 5 -, and have the general 
tight-binding form 

U aa ' 

Here bi >r7 (b\ a ) are destruction (creation) operators for 
a particle at site i in the internal state a, Jfj describes 
the coherent hopping of a particle from site i to site j 
(typically the nearest neighbor), and TJ°° describes the 
onsite (i = j) or offsite (i ^ j) two-body interactions 
between particles, with n i a = b\ a b ia . Hubbard models 
have a long history in condensed matter, where they 
have been used as tight-binding approximations of 
strongly correlated systems. For example, for particles 
being electrons in a crystal hopping from the orbital of 
a given atom to that of its nearest neighbor, tr repre- 
sents the electron spin. A (fermionic) Hubbard model 
comprising electrons in a 2D lattice with interspecies 
onsite interactions is thought to be responsible for 
the high-temperature superconductivity observed in 
cuprates 1 ^ 6 -. 

In recent years, Hubbard models have been shown to 
properly describe the low-energy physics of interacting 
bosonic and fermionic atoms trapped at the bottom of 
an optical lattic e 107 ' 108 . The resulting low-energy Hamil- 
tonians arc of the form of Eq. |(5J) above. Spectacular 
experiments with ultracold atoms have lead to the real- 
ization of the superffuid/Mott-insulator quantum phase 
transition for bosonic atom s 32 ' 33 , and great experimen- 
tal progresses with fermions promise to solve the phase 
diagram of the fermionic Hubbard model in 2D by per- 
forming an analog quantum simulation of Eq. ((SJ) with 
two-species cold fermion s 109 ' 110 . 

Since the interactions between cold atoms are short- 
ranged, in these systems Hubbard Hamiltonian typically 
have onsite interactions only [U^f in Eq. ([5])]. However, 



it has been shown that the presence of moderately 
long-range interactions in Eq. l[8]). such as nearest- 
neighbor interactions, can lead to interesting phases 
such as checkerboard solids and 2D supersolid a 87 ' 88 . 
Polar molecules in optical lattices can provide for offsite 
interaction s 82 ' 83 ' 84 ' 85 ' 86 which are strong, of the order of 
hundreds of kHz, and long-range, i.e. they decay with 
distance as l/\i — j| 3 . Due to these strong interactions, 
two molecules cannot hop onto the same site, and thus 
the particles are treated as effectively " hard-core" . 

An intriguing possibility is offered by the interaction 
engineering discussed above in the context of the real- 
ization of effective lattice models where particles interact 
via exotic (extended) Hubbard Hamiltonians. An exam- 
ple of this is given in Refi^, where it is shown how to 
engineer the following Hubbard-likc Hamiltonian 

h = -jj2btbj+j2^r n i n 3+ -j^r n i% n fc> ( 9 ) 

where Wijkriinjnk is an offsite three-body term. The 
latter is tunable independently of the two-body term 
UijUiUj, to the extent that it can be made to dominate 
the dynamics and determine the groundstate properties 
of the system. In contrast to the common approach to 
derive effective many-body terms from Hubbard models 
involving two-body interactions, which are obtained in 
a J <C U perturbation theory, and are thus necessarily 
smalUii, the derivation of the Hubbard model Eq. ^ 
is based directly on the effective many-particle potential 
Eq. ([2]). Thus, all the energy scales in Eq. ([9]) can be 
tuned independently, which allows to obtain compara- 
tively large hopping rates determining the time and tem- 
perature scales to observe exotic quantum phases. This 
is important, since e.g. in ID analytical calculations sug- 
gest that the Hamiltonian Eq. © has a rich groundstate 
phase-diagram, supporting valence-bond, charge-density- 
wave and supcrffuid phases 8 ^. In Sect. [TVl below we pro- 
vide the microscopic derivation of the effective interaction 
potentials of Eq. ^ . 

E. Lattice Spin models 

The Hamiltonian Eq. §5§ can be generalized to include 
other internal degrees of freedom for each molecule in 
addition to rotation. This offers new possibilities to en- 
gineer effective interactions and novel many-body phases. 
For example, the addition of a spin-1/2 (qubit) degree of 
freedom to polar molecules trapped into an optical lat- 
tice allows to construct a complete toolbox for the simula- 
tion of any permutation symmetric lattice spin models^ 6 -. 
Lattice spin models are ubiquitous in condensed matter 
physics where they are used as simplified models to de- 
scribe the characteristic behavior of more complicated 
interacting physical systems. 

The basic building block is a system of two polar 
molecules strongly trapped at given sites of an optical 
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FIG. 4: Example anisotropic spin models that can be simulated with polar molecules trapped in optical lattices, (a) Square 
lattice in 2D with nearest neighbor orientation dependent Ising interactions along x and z. Effective interactions between 
the spins Si and S2 of the molecules in their rovibrational ground states are generated with a microwave field E(t) inducing 
dipole-dipole interactions between the molecules with dipole moments Di and D2, respectively, (b) Two staggered triangular 
lattices with nearest neighbors oriented along orthogonal triads. The interactions depend on the orientation of the links with 
respect to the electric field. (Dashed lines are included for perspective.) 



lattice, where the spin-1/2 (or qubit) is represented by 
a single electron outside a closed shell of a 2 S 1 / 2 het- 
eronuclear molecule in its rotational ground state, as 
provided e.g. by alkaline-earth monohalogenides. As 
dicussed above, heteronuclear molecules have large per- 
manent electric dipole moments, which are responsible 
for strong, long-range and anisotropic dipole-dipole in- 
teractions, whose spatial dependence can be manipulated 
using microwave fields. Accounting for the spin-rotation 
splitting of molecular rotational levels these dipole-dipole 
interactions can be made spin- dependent. General lattice 
spin models are then readily built from these binary in- 
teractions. Although in this review we will present results 
for spin-1/2 models only, we notice that the inclusion of 
hypcrfine effects offers extensions to spin systems with 
larger spin. For example, the design of a large class of 
spin-1 interactions for polar molecules has been shown in 
Refi^, which allows e.g. for the realization of a general- 
ized Haldane model in 1D^. 

Two highly anisotropic models with spin-1/2 particles 
that can be simulated arc illustrated in Figs. [3Ja) and 
IDJb) respectively. The first takes place on a square 2D 
lattice with nearest neighbor interactions 

i-\ £-1 

<L=EE J KXi+i+ cos W+i,i)' ( 10 ) 
i=i 3=1 

Introduced by Doucot et al*^ in the context of Joseph- 
son junction arrays, this model (for £ 7^ ±7r/2) admits 
a 2- fold degenerate ground subspacc that is immune to 
local noise up to £-th order and hence is a good candidate 
for storing a protected qubit. 



The second, occurs on a bipartite lattice constructed 
with two 2D triangular lattices, one shifted and stacked 
on top of the other. The interactions are indicated by 
nearest neighbor links along the x,y and z directions in 
real space: 



h£L = j± E ^t+J± E °j°t+ J > E 

x— links y— links 2— links 

(11) 

This model has the same spin dependence and nearest 
neighbor graph as the model on a honeycomb lattice in- 
troduced by Kitae\iiii He has shown that by tuning 
the ratio of interaction strengths | J±\/\ J z \ on e can tune 
the system from a gapped phase carrying abclian anyonic 
excitations to a gapless phase which in the presence of a 
magnetic field becomes gapped with non-abelian excita- 
tions. In the regime | Jj_|/| J z \ -C 1 the Hamilonian can be 
mapped to a model with four body operators on a square 
lattice with ground states that encode topologically pro- 
tected quantum memorj*ii£. One proposalii^ describes 
how to use trapped atoms in spin dependent optical lat- 
tices to simulate the spin model H^ n . There the induced 
spin couplings are obtained via spin dependent collisions 
in second order tunneling processes. Larger coupling 
strengths as provided by polar molecules are desirable. 
In both spin models (I and II) above, the signs of the 
interactions are irrelevant although one is able to tune 
the signs if needed. 



8 



F. Hubbard models in self-assembled dipolar 
lattices 




FIG. 5: Floating lattices of dipoles: A self-assembled crystal 
of polar molecules with dipole moment d c provides a 2D pe- 
riodic honeycomb lattice V cp (darker shading corresponds to 
deeper potentials) for extra molecules with dipole d p <C d c 
giving rise to a lattice model with hopping J and long-range 
interactions Vij. 

In Hubbard models with cold atoms or molecules in 
optical lattices there is no phonon degrees of freedom 
corresponding to an intrinsic dynamics of the lattice, as 
the back action on the optical potentials is typically neg- 
ligible. Thus atomic and molecular Hubbard models al- 
low the study of strong correlations in the absence of 
phonon effects. However, the simulation of models where 
the presence of (crystal) phonons strongly affects the 
(Hubbard) dynamics of the particles remains a challenge. 
These models are of fundamental interest in condensed 
matter physics, where they describe polaronic and/or su- 
perconducting matcrialsiiS. In the context of atoms, one 
example is to immerse atoms moving on a lattice into 
a BEC of a second atomic species, representing a bath 
of Bogoliubov excitations. A second example is a self- 
assembled floating lattice of molecules as discussed in 
Sect. Ill Bt which provides a periodic potential for extra 
atoms or molecules, whose dynamics can be again de- 
scribed in terms of a Hubbard model^. Phonon degrees 
of freedom enter as vibrations of the dipolar lattice. 

The Hamiltonian for extra atoms or molecules in a self- 
assembled dipolar lattice is 

h = -j c ! c j + iE^ c H c j c ' 

<i,j> i,j 
1 1,3 

Here, the first and second terms define a Hubbard-likc 
Hamiltonian for the extra-particles of the form of Eq. (jSJ), 



where the operators Cj (cj) are destruction (creation) op- 
erators of the extra-particles. However, the third and 
fourth terms describe the acoustic phonons of the crys- 
tal and the coupling of the extra-particles to the crystal 
phonons, respectively. Here, a q destroys a phonon of 
quansimomentum q in the mode A. Tracing over the 
phonon degrees of freedom in a strong coupling limit 
provides effective Hubbard models for the extra-particles 
dressed by the crystal phonons 

H = -J ^ c\ Cj + \Y^ Vijcl^CjCi 

<i,j> i,j 

The hopping of a dressed extra-particle between the min- 
ima of the periodic potential occurs at a rate J, which 
is exponentially suppressed due to the co-propagation of 
the lattice distortion, while offsite particle-particle inter- 
actions Vij are now a combination of direct particle- 
particle interactions and interactions mediated by the 
coupling to phonons. The setup we have in mind is 
depicted in Fig. [^b), where extra-particles which are 
molecules with a dipole moment d p <C d c interact repul- 
sively with the crystal molecules, and thus see a periodic 
(honeycomb) lattice potential. 

The distinguishing features of this realization of lattice 
models arc: (i) Dipolar molecular crystals constitute an 
array of microtraps with its own quantum dynamics rep- 
resented by phonons (lattice vibrations), while the lattice 
spacings are tunable with external control fields, ranging 
from a /im down to the hundred nm regime, i.e. poten- 
tially smaller than for optical lattices, (ii) The motion 
of the extra particles is governed by an interplay of Hub- 
bard (correlation) dynamics in the lattice and coupling 
to phonons. The tunability of the lattice allows to access 
a wide range of Hubbard parameters and phonon cou- 
plings. Compared with optical lattices, for example, a 
small scale lattice yields significantly enhanced hopping 
amplitudes, which set the relevant energy scale for the 
Hubbard model, and thus also the temperature require- 
ments for realizing strongly correlated quantum phases. 



III. ENGINEERING OF INTERACTION 
POTENTIALS 

In this section we show in some detail how to realize 
a collisionally stable 2D setup where particles interact 
via a purely repulsive 1/r 3 potential, by using a combi- 
nation of a DC field and of tight optical confinement in 
the field's direction. We then sketch how to design more 
complicated interactions using a combination of AC, DC, 
and optical fields, by focussing on the step-like potential 
of Fig.QJb). This engineering of interaction potentials is 
at the core of the realization of the strongly correlated 
phases and quantum simulations discussed in Sect. IfVI 
below. 
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A. Molecular Hamiltonian 

We consider spin-less polar molecules in their elec- 
tronic and vibrational ground-state, with spectral nota- 
tion X 1 S(0). In the following, we are interested in ma- 
nipulating their rotational states using DC and AC elec- 
tric fields and in confining their motion using a (optical) 
far-off-resonance trap (FORT). The application of these 
external fields will serve as a key element to engineer ef- 
fective interaction potentials between the molecules. 

The low energy effective Hamiltonian for the exter- 
nal motion and internal rotational excitations of a single 
molecule is 

P 2 

H{t) = 7^ + #rot + H DC + H AC (t) + ffopt(r), (12) 

where p 2 /2m is the kinetic energy for the center-of-mass 
motion of a molecule of mass m, H mt accounts for 
the rotational degrees of freedom, while the terms 
Hoc, HAc{t) an d i/ op t(r) refer to the interaction 
with electric DC and AC (microwave) fields and to 
the optical trapping of the molecule in the ground 
electronic-vibrational manifold, respectively. In the 
following we consider tight harmonic optical traps with 
a frequency lo±_ = 2n x 150kHz, which is the same 
for all the relevant rotational states of the molecule. 
That is, we neglect possible tensor-shifts induced by 
the optical trapping in the energies of the excited 
rotational states of the molecules, which in general 
can be compensated for by an appropriate choice of 
additional laser fields^. Thus, for a confinement along 
e z , if opt (r) reads H opt (r) = muj 2 L z 2 /2, independent of 
the internal (rotational) 44 ' 119 . 

Rotational spectrum :- The term H rot in Eq. (| 12[) is the 
Hamiltonian for a rigid spherical rotopi 2 ^ 

H mt =BJ 2 , (13) 

which accounts for the rotation of the internuclear axis 
of a molecule with total internal angular momentum 
ji20 1 i2i 1 i22 Rotations are the lowest-energy internal 
excitations of the molecule. Here B is the rotational 
constant for the electronic-vibrational ground state, 
which is of the order of B ~ h 10 GHai 2 ^. We denote 
the energy eigenstates of Eq. (fP3"j) by | J, M) , where J is 
the quantum number associated with the total internal 
angular momentum and M is the quantum number asso- 
ciated with its projection onto a space-fixed quantization 
axis. The excitation spectrum is Ej = BJ(J+ 1), which 
is anharmonic. Each J-level is (2 J + l)-fold degenerate. 

A polar molecule has an electric dipole moment, d, 
which couples its internal rotational levels and for S- 
molecules is directed along the internuclear axis e a b, i.e. 
d = de a b. Here, d is the "permanent" dipole moment 
of a molecule in its electronic-vibrational ground-state. 
This dipole moment is responsible for the dipole-dipole 
interaction between two molecules. 



The spherical component d q = e q ■ d of the dipole oper- 
ator on the space-fixed spherical basis {e_i,eo,ei}, with 
e g= o = e z a nd e±i = ^(e-x ± ie y )/>/2 couples the rota- 
tional states | J, M) and | J ± 1, AI + q) according to 

(J ± 1, M + q\d q \J, M) = d(J, M; 1, q\ J ± 1, M + q) x 

x(J,0;l,0|J ± l,0)^±L:, 

where ( J\, M\\ J 2 , Af 2 | J, M) arc the Clebsch-Gordan- 
cocfficicnts. This means that for a spherically-symmetric 
system the eigenstates of the rotor have no net dipole- 
moment, (J, M\d\J, M) = 0. However, the dipole 
coupling to an external electric field breaks this spherical 
symmetry by aligning each molecule along the field's 
direction. This induces a dressing of the rotational 
energy levels of the molecule, and a corresponding finite 
dipole moment in each rotational state, as explained 
below. 

Coupling to external electric fields :- The terms Hy>c 
and BiAc(t) in Eq. (fT2"]) are the electric dipole interac- 
tion of a molecule with an external DC electric field 
Edc = £i>c e z directed along e = e z , and with AC 
microwave fields Eac(*) = EAce~" luJt e q + c.c, which are 
linearly (q = 0) or circularly polarized (q = ±1) relative 
to e 2 , respectively. Here we have neglected the spatial 
dependence of Eac since in the following we are inter- 
ested in dressing the rotational states of the molecules 
with microwave fields, whose wavelengths are of the or- 
der of centimeters, and thus much larger than the size of 
our system. Then, the terms i?DC and Hac (t) in Eq. (fT2|) 
read 

H DC = -d-E DC = -d E DC , (14a) 
H AC (t) = -d-E AC (t) = -dqEAce-™* + h.c. 

(14b) 



In the presence of a single DC electric field Edc 
(EaCj w± = 0), the internal Hamiltonian is that of a rigid 
spherical pendulum^ H = H rot + Hue = BJ 2 — do^DC , 
which conserves the projection of the angular momen- 
tum J on the quantization axis, i.e. M is a good quan- 
tum number. Thus, the energy eigenvalues and eigen- 
states are labeled as Ej.m and \4>j.m), respectively, where 
each eigenstate \4>j,m) is a superposition of various states 
I J, M) mixed by the electric dipole interaction. 

The effects of a DC electric field on a single polar 
molecule are shown in figure Fig. and they amount 
to: (a) split the (2J + l)-fold degeneracy in the rotor 
spectrum, and (b) align the molecule along the direction 
of the field. The latter corresponds to inducing a finite 
dipole moment in each rotational state. For weak fields 
(3 = dE^c/B <C 1, the state \4>j,m) and its associated 
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FIG. 6: Solid lines: Energies Ej t M (left) and states \4>j,m} 
(right) with J — 0, 1, for a molecule in a weak DC electric 
field Edc = -B D ce with (3 = dE D c/B < 1. The DC-field- 
induced splitting fi<5 and the average energy separation to are 
hS = 3B/3 2 /20 and to = 273 + B/? 2 /6, respectively. Dashed 
and dotted lines: Energy levels for a molecule in combined DC 
and AC fields (The AC-Stark shifts of the dressed states are 
not shown). Dashed line: The AC field is monochromatic, 
with frequency uo, linear polarization q = 0, and detuning 
A = uj — (ui + 25/3) > 0. Dotted lines: Schematics of energy 
levels for an AC-field with polarization q = ±1 and frequency 
uj' 7^ UJ. 



induced dipole moment approximately read 



>J,M 



2^3(27+1)' 



p V(J+1) 2 -A/ 2 
"2^(J+1)3( 2 J+1) 



\J+l,M), (15) 



and 



>J,M 



d|0j, 



d/3 



3M 2 /J(J+ 1) - 1 
(2J- l)(2J + 3) 



e , 



respectively. Thus, the ground state acquires a finite 
dipole moment (0o,oMo|0o,o) = d/3/3 along the field axis, 
which is at the origin of ground-state dipole-dipole inter- 
actions between polar molecules. 

We notice that the for a typical rotational constant, 
B ~ h 10 GHz, and a dipole-moment d ~ 9 Debye the 
condition /3 <C 1 corresponds to considering DC fields 
(much) weaker than B/d ~ 2 kV/cm. 

Individual transitions of the internal Hamiltonian can 
be addressed by applying one (or several non-interfering) 
microwave field Eac(£)i which can be linearly or cir- 
cularly polarized. This is shown in Fig. [H] for transi- 
tions coupling the J = and J = 1 manifolds, where 
the Rabi frequency and the detuning A are Q = 
EAc(<j>i,q\d q \<po,o}/K and A = uj - (E 1<q - E . )/h, re- 
spectively 



Dressed energy levels of a molecule are obtained by di- 
agonalizing the Hamiltonian H = H Iot + i?DC + ^ac (t) 
in a Floquet picture. That is, first, the Hamilto- 
nian is expanded on the basis \4>j,m), which diagonal- 
izes the time- independent part of H as H ro t + Hdq = 
S 7 a/ \4 > J,m)Ej^m(4 ) J,m\i and then the time-dependent 
wave-function is expanded in a Fourier series in the AC 
frequency uo. After applying a rotating wave approxi- 
mation, i.e. keeping only the energy conserving terms, 
one obtains a time-independent Hamiltonian H, whose 
eigenvalues correspond to the dressed energy levels^. 



B. Two molecules 

We now consider the interactions of two polar 
molecules j = 1,2 confined to the x — y plane by a tight 
harmonic trapping potential of frequency directed 
along z. The interaction of the two molecules at a dis- 
tance r = Y2 — t\ = re r is described by the Hamiltonian 



(16) 



where Hj (t) is the single-molecule Hamiltonian Eq. (fT! 
and Vdd(r) is the dipole-dipole interaction of Eq. (H|). 

In the absence of external fields £bc = -^ac = 0, 
the interaction of the two molecules in their rotational 
ground state is determined by the van-der-Waals attrac- 
tion Vvdw ~ C^/r 6 with Cq rts —d 4 /6B. This expression 
for the interaction potential is valid outside of the 
molecular core region r > rs = (d 2 /B) 1 ^ 3 , where rs 
defines the characteristic length where the dipole-dipole 
interaction becomes comparable to the splittings of the 
rotational levels. In the following we show that it is 
possible to induce and design interaction potentials 
which are long-range, by dressing the interactions with 
appropriately chosen static and/or microwave fields. In 
fact, the combination of the latter with low-dimensional 
trapping allows to engineer effective potentials whose 
strength and shape can be both tuned. The derivation of 
the effective interactions proceeds in two steps: (i) We 
derive a set of Born-Oppcnheimer (BO) potentials by 
first separating Eq. (|16p into center-of-mass and relative 
coordinates, and diagonalizing the Hamiltonian iJ re i 
for the relative motion for fixed molecular positions. 
Within an adiabatic approximation, the corresponding 
eigenvalues play the role of an effective 3D interac- 
tion potential in a given state manifold dressed by 
the external field. (ii) We eliminate the motional 
degrees of freedom in the tightly confined z direction to 
obtain an effective 2D dynamics with interaction V^g(p). 

In the following we show how to design interaction po- 
tentials, presenting in some details the simplest case of 
a purely repulsive 1/r 3 potential in 2D, obtained using 
a static electric field (Sect. IHI B ip . We then sketch how 
to design more elaborate potentials using a combination 
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of static and microwave fields coupling the lowest rotor 
states of each molecule (Sect. MI B 21) . 



1. Designing the repulsive 1/r 3 potential in 2D 

Collisions in a DC field:- We consider a weak static 
electric field applied in the ^-direction E = i?Dc e o with 
[3 = dEnc/B <C 1, and in the absence of optical trap- 
ping = 0). The effective interaction potentials for 
the collision of the two particles can be obtained in the 
adiabatic approximation by neglecting the kinetic energy 
and by diagonalizing the following Hamiltonian H re \ for 
fixed particle positions 

2 

H iel = i Bj2 j - E ncd 0;j ] + V dd (r) 



3 = 1 



£|* n (r)>25 B (r)<* B (r)|, 



(17) 



where E„(r) and |$„(r)) arc the nth adiabatic energy 
eigenvalues and two-particle eigenfunctions, respectively. 
In the limit r —* oo the latter are symmetrized products 
of the single-particle states l^j^M^j of Eq. (|15p . while for 
finite r they are superpositions of several single-particle 
states, which are mixed by the dipole-dipole interaction 
Vdd(r). The quantity n = (J;M;o~) is the collective 
quantum number labeling the eigenvalues E n (r), with 
J = Ji + J2 the total number of rotational excitations 
shared by the two molecules, M = \M\\ + IM2I the total 
projection of angular momentum onto the electric field 
direction, and a = ± the permutation symmetry asso- 
ciated with the exchange of the two particles. We note 
that, because of the presence of the DC field, here J is a 
simple label for the various energy manifolds, and not a 
quantum number. 

Since we are mainly interested in ground-state colli- 
sions, in the following we restrict our discussion to the 
Jj = and 1 manifolds of each molecule, which amounts 
to taking into account 16 rotational two-particle states. 
Figure [7] shows the corresponding eigenvalues E n (r) as a 
function of the interparticle distance r, for (3 = 1/5. The 
vector r is expressed in spherical coordinates r = (r, tp), 
with ■& and <p the polar and azimuthal angles, respec- 
tively, and z = r cos fl. Figure Eta) shows that the en- 
ergy spectrum has a markedly different behavior in the 
molecular core region r < tb and for r > rg, with 
tb = (d 2 /B) 1 / 3 . In fact, for r <tb the energy spectrum 
is characterized by a series of level crossings and anti- 
crossings, which make the satisfaction of the adiabatic 
approximation generally impossible. For r > tb the en- 
ergy levels group into well-defined manifolds, which are 
approximately spaced by an energy 2B, corresponding to 
a quantum of rotational excitation. In the following we 
focus on this region r > rs, where the adiabatic approx- 
imation can be fulfilled. 

Figures [7](b,c) and Figs.[7fd,e) are blow-ups of the two 
lowest-energy manifolds of Fig.[7]Ja) in the region r > tb, 



for t3 = tt/2 and # = 0, respectively. Figure. [7jb) and 
Fig. EJd) show that the excited state manifold with one 
quantum of rotation ( J\ + 3% = 1) is asymptotically split 
into two sub-manifolds. This separation corresponds to 
the electric-field-induced splitting of the Jj = 1 manifold 
of each molecule, and it is thus given by HS = 3B0 2 /20, 
see caption of Fig. [51 More importantly, Figs.[7Jc) and (e) 
show that the effective ground-state potential Eq (r) has a 
very different character for the cases ■& — tt/2 and i9 = 0, 
respectively. In fact, in the case $ = n/2 [Fig. E£c)], 
corresponding to collisions in the (z = 0)-plane, the po- 
tential is attractive for r < r*, while for r > it turns 
into repulsive and it decays at large distances as 1/r 3 , 
where is a characteristic length to be defined below. 
On the other hand, for 1? = [see Fig. EJe)] the potential 
is purely attractive, with dipolar character. This change 
in character of the ground-state potential as a function 
of $ is captured by the following analytic expression for 
-Eo ; o;+( r ), as derived in perturbation theory in Vdd(r)/B, 



V c f(r) ee £ ;0;+(r) 



C3 

72 al , 



1 - 3 cos 2 ■&) + 



(18) 



Here, the constants C3 ~ d f3 /9 and Cq w —d /QB are 
the dipolar and van-der-Waals coefficients for the ground- 
state, respectively, and the constant term 2£'o,o = 
—f3 2 B/3 due single-particle DC Stark-shifts has been 
neglected. Equation (fl8)) is valid for r 3> tb and 



Vdd(r)/£> <C 1, and it shows that the potential V^P(r) 



has a local maximum in the plane z 
position r*, defined as 



r cos i3 — at the 



2\C e 



1/3 



/3d 2 



\B[3' 2 



1/3 



(19) 



where the dipolc-dipolc and van-der-Waals interactions 
become comparable. The height of this maximum is 



c 3 2 

4|C 6 | 



B(3 4 
~54~ 



(20) 



and the curvature of the potential along z [d 2 V(r 



r*,z = 0) = -6C 3 /r 5 , = 
istic frequency 



?/2] defines a character- 



•jJr 



1/2 



(21) 



to be used below. The latter has a strong dependence 
/3 8 / 3 = (dE-oc/B) 8 / 3 on the applied electric field. 
For distances r r± the dipole-dipole interaction 
dominates over the van-der-Waals attractive potential, 
and V c 3 ff D (r) - C 3 (l - 3cos 2 tf)/r 3 , see Ref^. Thus, 
if it were possible to confine the collisional dynamics 
to the (z = 0)-plane with r 3> r^,rs, purely repulsive 
long-range interactions with a characteristic dipolar 
spatial dependence ~ 1/r 3 could be attained. In the 
following we analyze the conditions for realizing this 
setup, using a strong confinement along z as provided, 
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FIG. 7: BO-potentials Ej-M-,a(r, •&) for two molecules colliding in the presence of a DC field, with f3 = dEuc/B = 1/5. The solid 
and dashed curves correspond to symmetric (a — +) and antisymmetric (a — — ) eigenstates, respectively, (a): BO-potentials 
for the 16 lowest-energy eigenstates Ej- t M\<r{r, The molecular-core region is identified as the region r < tb — (d 2 / B) 1 ^ 3 , 
while for r ^> tb the eigestates group into manifolds separated by one quantum of rotational excitation 2B. (b) and (d): 
Blow-ups of the first-excited energy manifold of panel (a) in the region r > tb for ■& = n/2 and ■£) — 0, respectively. Note 
the electric-field-induced splitting HS = SB f3 2 /20. The distance r$ where the dipole-dipole interaction becomes comparable to 
HS is rg = (d 2 /hS) 1 ^ 3 . (c) and (e): Blow-ups of the ground-state potential Eo > o-+(r, ■&) of panel (a) in the region r > t\b for 
■& — tt/2 and ■& = 0, respectively. The distance r* of Eq. (|19|l . where the dipole-dipole interaction becomes comparable to the 
van-der-Waals attraction is indicated. Note the repulsive (attractive) character of the potential for # = tv/2 = 0) and r > r*. 



e.g. by an optical trapping potential. 

Parabolic confinement :- The presence of a finite trap- 
ping potential of frequency lo±_ in the z-direction provides 
for a position-dependent energy shift of Eq. (jTHJ) . The 
new potential reads 

V(r) = ^(l-3co S ^) + ^+ 1 -m^. (22) 

As noted before, for z = the repulsive dipole-dipole 
interaction dominates over the attractive van-der-Waals 
at distances r 3> r* given in Eq. (fTi?)) . In addition, for 
u>± > the harmonic potential confines the particle's 
motion in the z direction. Thus, the combination of the 
dipole-dipole interaction and of the harmonic confine- 
ment yields a repulsive potential which provides for a 
three-dimensional barrier separating the long-distance 
from the short-distance regime. If the collisional energy 
is much smaller than this barrier, the particle's motion is 
confined to the long-distance region, where the potential 
is purely repulsive. 

Figure [8] is a contour plot of V(r) in units of V*, for 
(3 > and ui± = w c /10, with r = (p, z) = r(sin$,cosi9) 
(the angle ip is neglected due to the cylindrical sym- 
metry of the problem). Darker regions correspond to 



a stronger repulsive potential, while the white region 
at p f» is the short-range, attractive part of the 
potential. The repulsion due to the dipole-dipole and 
harmonic potentials is distinguishable at ~ 

and 7, respectively. The lesser-dark regions located 
at (p±,±z±) = ^j_(sin^j_, ± cosi9j_) correspond to the 
existence of two saddle points positioned in between 
the maxima of V(r), with £±_ = (^C^/mt^) 1 / 5 and 
cos^j^ = y/l — (r±/£±) 3 /V5, [sec circles in Fig. [8]. 
These saddle points act as en effective potential barrier 
separating the attractive part of the potential present at 
r < lj_ from the region r 3> £±_ > r*,rs where the ef- 
fective interaction potential Eq. (|22[) is purely repulsive. 
For collisional energies smaller than the height of this 
barrier the dynamics of the particles can be reduced to 
a quasi two-dimensional (2D) dynamics, by tracing over 
the fast particle motion in the z-direction. We notice 
that the existence of two saddle points at distances 
r ~ £±_ separating the long- from the short-distance 
regimes is a general feature of systems with a compar- 
atively weak transverse trapping uj±/lu c < 1, with u> c 
defined in Eq. (|21j) . In fact, for a strong transverse 
trapping lu±>ui c the two saddle points collapse into 
a single one located at z = 0, and p — £± ~ r*. In 
this limit the dynamics is purely 2D, with the particles 
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FIG. 8: (a) Contour plot of the effective potential V(p,z) of Eq. (|22[) . for two polar molecules interacting in the presence 
of a DC field j3 > 0, and a confining harmonic potential in the z-direction, with trapping frequency u>± — w c /10, where 
lu c = {12C 3 /mr r l) 1/2 of Eq. J2lj and r* = (2|C 6 1/C*3) 1/3 of Eq. JT9J. The contour lines are shown for V(p,z)/V* > 0, with 
V* = i?/9 4 /54. Darker regions represent stronger repulsive interactions. The combination of the dipole-dipole interactions 
induced by the DC field and of the harmonic confinement leads to realizing a 3D repulsive potential. The repulsion due to 
the dipole-dipole interaction and of the harmonic confinement is distinguishable at z ~ and a/r* ~ ±7, respectively. Two 
saddle points (circles) located at (p±,±z±) separate the long-distance region where the potential is repulsive ~ l/r' ! from 
the attractive short-distance region. The gradients of the potential are indicated by dash-dotted lines. The thick dashed line 
indicates the instanton solution for the tunneling through the potential barrier, (b) The euclidian action Se as a function of 
uj±/u) c (solid line). For oj±_ < uj' c ~ 0.88 uj c (wi > uj' c ) the "bounce" occurs for z(0) 7^ (within the plane z(0) = 0), see text. 
The point u>' c is signaled by a circle. For uj± > ui' c the action is 5*e ~ 5.785*o, with So = ^Jm\C§\ / '%r\ , which is lu ^-independent, 
consistent with the "bounce" occurring in the (z — 0)-plane (see text). 



strictly confined to the (z = 0)-plane. 

Collisional stability :- When an ensemble of polar 
molecules is considered, inelastic collisions and three 
body recombination may lead the system to a potential 
instability, associated with the attractive character of 
the dipole-dipole intcra<:tionSLfi&fimZUmZiJ&. In 
our discussion, this instability is associated with the 
population of the short-distance region r < £±, which 
can be efficiently suppressed for strong dipole-dipole 
interactions and transverse confinement. In fact, for 
collisional energies smaller than the potential barrier 
V(p±,±z±) the particles are mostly confined to the 
long-distance regime, where they scatter elastically. 
That is, when a cold ensemble of molecules is considered 
the barrier provides for the stability of the system by 
"shielding" the short-distance attractive part of the 
two-body potential. In this limit, residual losses are due 
to the tunneling through the potential barrier at a rate 
r, which can be efficiently suppressed for reasonable 
values of (3 and u>± , as shown below. 

The tunneling rate T = T e~ SE ^ h through the 
barrier V(pj_, ±zj_) can be calculated using a semi- 
classical/instanton approacbA2i. The euclidian action 
Se, which is responsible for the exponential suppression 



of the tunneling, is plotte d in F ig. [5] as a function of 
uj_i_/uj c , in units of So = \Jm\C§\/r1. The figure shows 
that Se has different behaviors for luj_ -C w c and luj_ 3> 
lu c . In fact, for uj±_ <C lu c the action increases with in- 
creasing cj_l as S E « 7.015" (w_ l /cj c ) 1 / 5 = 1.43ft(£ ± /a ± ) 2 
(dotted line), which depends on the confinement along z, 
via a±_ = {h/muj^) 1 / 2 . On the other hand, for uj±_ uj c 
it reads ~ 5.78 So, which is ^-independent. The 
transition between the two different regimes mirrors the 
change in the nature of the underlying potential V(r). 
In particular, for u±_ > lu c the dynamics is strictly con- 
fined to the plane z = and thus it becomes indepen- 
dent of luj_. The constant Tq is related to the quan- 
tum fluctuations around the semiclassical trajectory, and 
its value is strongly system-dependent. For the crys- 
talline phase of Rcf.— , it is the collisional " attempt fre- 
quency" , proportional to the characteristic phonon fre- 
quency Tq ~ y^Cs/ma 5 , with a the mean interparticle 
distance. 

In the limit of strong interactions and tight transverse 
confinement T rapidly tends to zero. We illustrate 
this for the example of SrO, which has a permanent 
dipole-moment of d ks 8.9 Debye and mass m = 104 amu. 
Then, for a tight transverse optical lattice with harmonic 
oscillator frequency uj± = 2ir x 150kHz and for a DC-field 
= dE DC /B = 1/3 we have (C*|m 3 w J ./8S 5 ) 1 /5 ~ 3.39 
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and wc obtain T/r ps e - 5 - 86x3 - 39 ps 2 x 10" 9 . Even 
for a DC field as weak as /3 = 1/6 we still ob- 
tain a suppression by five order of magnitudes, as 
T/r « e -5.86xi.94 _ 1Q -5 Thig calculation confirms 

that a collisionally stable setup for polar molecules 
in the strongly interacting regime can be realized by 
combining the strong dipole-dipole interactions with a 
tight transverse confinement. 

Effective 2D interaction:- The effective two- 
dimensional interaction potential is obtained by 
integrating out the fast particle motion in the trans- 
verse direction z. For r > £± 3> a±, the two-particle 
eigenfunctions in the z-direction approximately factorize 
into products of single-particle harmonic oscillator wave- 
functions tpk 1 (zi)tpk 2 ( z 2), and in first order perturbation 
theory in V 2 ® /hu± the effective 2D interaction potential 
V 2 ® reads 

V™(p) « - JU J dze-*'**.Vg{v). (23) 
For large separations p S> t± the 2D potential reduces to 

t/2D / \ ^3 

Kg (p) = — . 

which is a purely repulsive 2D interaction potential. The 
derivation of V 2 q(p) is one of the central results of this 
section. We show below (Sect. IIV[) that the use of this 
interaction potential leads to the realization of interesting 
many-body phases, in the context of condensed matter 
applications using cold molecular quantum gases. 

2. Designing ad-hoc potentials with AC-fields 

Above we have shown how to design 2D effective 
groundstate interactions which are purely repulsive and 
decay as ~ 1/r 3 . The use of one or several non interfering 
AC fields allows to engineer more complicated interac- 
tions by combining the spatial texture of the adiabatic 
groundstate potential of the two-particle spectrum with 
that of selected excited potentials, in a dressed picture. 
This mixing of ground and excited-state potentials is 
favored by the dipole-dipole interactions which split 
the degeneracy of the excited-state manifolds of the 
two-particle spectrum and render the state-selectivity of 
the AC fields space- dependent, as explained below. In 
combination with a strong optical confinement, and due 
to the long lifetimes of the excited rotational states^, 
this allows for the realization of collisionally stable 
setups for molecules in the strongly interacting regime. 

We exemplify the situation above by considering the 
case of a single AC field Eac(*) = Ex C e~ lut e q + c.c. 
which is added to the configuration of Fig. [7] (interactions 
in the presence of a static electric field Edc = PBe z ). 
The field's polarization is chosen to be linear (q = 0) 



and the frequency uj is blue-detuned from the (|0o,o) ~~ * 
|(/)i j o))-transition of the single-particle spectrum by an 
amount A = uj - 2B/H > 0. The effects of the AC- 
field on the two-particle scattering can be summarized as: 
(a) Inducing oscillating dipole-moments in each molecule, 
which determine long-range dipole-dipole interactions [in 
addition to those determined by the static field EdcL 
whose sign and angular dependence are given by the po- 
larization q; (b) Inducing a coupling of the ground and 
excited state manifolds of the two-particle spectrum at a 
resonant (Condon) point rc = (d 2 /3ft- A) 1 / 3 , where the 
dipole-dipole interaction becomes comparable to the de- 
tuning A. This coupling is responsible for an avoided 
crossing of the field-dressed energy levels at rc, whose 
properties depend crucially on the polarization q. This 
fact is at the core of the engineering of interaction poten- 
tials, in that the 3D effective dressed adiabatic ground- 
state interaction potential inherits the character of the 
bare ground and excited potentials for r 3> rc and 
r <C rc, respectively. 

The setup above is illustrated in Fig. [D(a) and (b), 
where the continuous and dashed lines arc the bare 
(Eac = 0) symmetric and anti-symmetric potentials 
Ej;M;<,{y) of Fig. [3 respectively, and the presence of the 
AC-ficld is signaled by a black arrow at the resonant Con- 
don point rc- The presence of the weak DC field splits 
asymptotically the (J = l)-manifold by an amount HS as 
in Fig. [7Kb), allowing for a simple fulfillment of the adi- 
abatic approximation in the excited-state manifold for 
distances r 3> r$ = (d 2 /M) 1 / 3 . In fact, the energy of 
the £Ji ; o ; +(r) potential becomes degenerate with the en- 
ergy of other bare symmetric potentials only at distances 
r <C rg. In addition, we notice that the presence of the 
splitting hS also shifts the level crossing with antisym- 
metric states to small distances r<fj. 

For distances r 3> r$ = (d 2 /hS) 1 ^ 3 , we are allowed to 
consider only the four states of Fig. [D(b), since all the 
other potentials of the (J = 1)- and (J = 2)-manifolds 
are far detuned by an amount which is (at least) of 
order S » A and they are not coupled by the AC-field 
to the bare ground state £t),0;+( r )> due to the choice of 
field's polarization. Figure [^b) shows that the splitting 
induced by the dipole-dipole interaction in the (J = 1)- 
manifold renders the detuning A position-dependent, so 
that at rc the energy of the bare ground-state and that 
of the symmetric bare excited state become degenerate. 
The resulting dressed ground-state potential is sketched 
in Fig. [UJb) (thick black line) and it roughly corresponds 
to the bare E 0t0[+ (r) and Ei :0 -,+ (r) potentials for r > rc 
and r < rc, respectively. Accordingly, Fig. Etc) shows 
that the dressed groundstate potential Eo : Q. + (r), which 
has the highest energy, turns from weakly to strongly 
repulsive for r > rc and r rc, respectively. This 
change in the character of the ground-state interaction 
potential corresponds to the design of a "step-like" in- 
teraction. This example shows that the 3D ground-state 
interaction for two molecules can be strongly modified 
by the combined use of AC and DC fields, which is 



15 



[Ej, M , rT {r^/2)-2E Qfi ]/h 

^2:0:- 




[Ej. M , a (r,n/2)-2E 0fi ]/hA 



10 20 30 40 50 
rj r/r B r c 



10 20 30 40 50 60 70 80 90 100 o 
rs r c t/tb 







-C/0;0; + 










r 






10 20 30 40 50 60 70 80 90 100 
r ° r/r B 



FIG. 9: (a) Schematic representation of the effects of a DC and an AC microwave fields on the interaction of two molecules. 
The solid and dashed lines are the bare potentials E n (r) = Ej- t M;a(r,'d) of Sec. MI B II with i9 = n/2 for interactions in the 
presence of the DC field only, for the symmetric (a = +) and antisymmetric (a — — ) states, respectively. The DC field induces 
a splitting hS of the first-excited manifold of the two-particle spectrum. A microwave-field of frequency uj — aJ + 25/3 + Ais blue 
detuned by A > from the single-particle rotational resonance. The dipole-dipole interaction further splits the excited-state 
manifold, making the detuning space-dependent. Eventually, the combined energy of the bare ground-state potential i?o;0;+( r ) 
and of an AC photon (black arrow) becomes degenerate with the energy of the bare symmetric £7i ; o ; +(r,7r/2). The resonant 
point rc = (d 2 /3fiA) 1//3 occurs at r « 46 tb- (b) Blow-up of the potentials of panel (a) with M = 0. The dressed ground- 
state potential is sketched by a thick solid line, (c) The four potentials of panel (b) in the field-dressed picture. The dressed 
ground-state potential Eo;0;+(r,n/2) has the largest energy and is indicated by a thick solid line. 



the central result of this section. More complicated 
potentials can be engineered using multiple AC fields 
and different polarizations. 

Analogous to the case (Eac = 0) of Sect. 1MB II the 
interaction potential of Fig. [9] is actually repulsive along 
certain directions (e.g. 8 = tt/2, as shown in the figure), 
while it turns into attractive along others (e.g. 6 = 0, not 
shown). As for the (Eac = 0)-case of Sect. MI B II when 
more than two particles are considered this attraction 
can lead to many-body instabilities. Moreover, here the 
dressed potential i?o ; o;+( r ) of Fig. [91c) is not the lowest- 
energy potential, which in general can introduce addi- 
tional loss channels. The latter correspond to diabatic 
couplings to symmetric states for particles approaching 
distances r < rc, and are therefore present even in 
the simple two-particle collisional process, and to cou- 
plings to anti-symmetric states, which can be induced e.g. 
by three-body collisions or by non-compensated tensor- 
shifts for two optically-trapped particles. The presence 
of all of these loss channels may render impractical the 
realization of collisionally stable setups for strongly in- 
teracting molecular gases (although Refi^ for a solution 
involving the use of a circularly polarized AC field). How- 
ever, we have seen above that for the setup of Fig. [9] the 
presence of the static field shifts the various resonance 
points with the potentials which are responsible for these 
loss channels in the ($ = n/2)-plane (z = 0) to distances 



r -C rc- This suggests that by confining the particles 
motion to the plane z = by using a strong optical trans- 
verse confinement analogous to that of Sect. IIII B II it is 
possible to realize collisionally stable setups in the region 
r > rc ■ This scheme has been shown to work in Refi^ 
and thus the main message here is that a judicious combi- 
nation of the dipole-dipole interactions and of the optical 
confinement can act as an effective " shield" of the region 
r < rc where losses occur and thus the collisional setup 
can be made stable. The use of the step-like potential 
above and of other engineered potentials can lead to the 
realization of interesting phases for an ensemble of polar 
molecules in the strongly interacting rcgim o 50 ' 66 . 

IV. MANY-BODY PHYSICS WITH COLD 
POLAR MOLECULES 

A. 2D Self-Assembled crystals 

The above discussion of the intermolecular potentials 
and of the stability of collisional setups in reduced dimen- 
sionality provides the microscopic justification for study- 
ing an ensemble of polar molecules in 2D interacting via 
(modified) dipole-dipole potentials. At low temperatures 
T < the general many body Hamiltonian has the 

form of Eq. As an example of the possibilities of- 
fered by potential engineering to realize novel many-body 
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FIG. 10: (a) Quantum phases of 2D dipoles: Contour plot of the interaction strength ra = Dm/h 2 a as a function of the dipole 
moment d (in Debye) and of the interparticle distance a (in /im), with m the mass of a molecule (in atomic units 200xamu). 
The regions of stability of the 2D superfluid and crystalline phases where Huj± > D/a 3 are indicated, with ui±_ — 2ty x 150kHz 
the frequency of the transverse confinement, (b) PIMC-snapshot of the mean particle positions in the crystalline phase for 
N = 36 at rd ~ 26.5. (c) Density-density (angle-averaged) correlation function g2(r), for TV = 36 at Td ~ 11.8. (d) Superfluid 
density p 3 and (e) static structure factor S(K)/N as a function of r<j, for N — 36 (circles) and N = 90 (squares). 



quantum phases, we here focus on bosonic particles in- 
teracting via the effective potential 

V™(p) = D/p 3 (24) 

derived in Sect. MI B 11 This is the simplest attainable in- 
teraction potential, and thus the one which is experimen- 
tally most interesting in the short term. However, despite 
the simplicity of the interaction potential, the Hamil- 
tonian Eq. ([SJ gives rise to novel quantum phenomena, 
which have not been accessed so far in the context of cold 
neutral atoms and molecules. In particular, by means of 
Path Integral Monte-Carlo simulations (PIMC), in Rcf. 50 
we show the appearance of a self-assembled crystalline 
phase and an associated quantum melting transition into 
a superfluid as a function of the interaction strength r c i- 
As explained in Sect. [HI the latter is the ratio 



between the interaction energy D/a 3 and the kinetic 
energy h 2 /ma 2 at the mean interparticle distance a, 
with m the mass of a molecule. 

In Fig. a tentative phase diagram is sketched for 
the two-dimensional system of bosonic dipoles. In the 
limit of weak interactions r<j < 1, the ground state is a 
superfluid (SF) with a finite (quasi) condensate. The SF 
is characterized by a superfluid fraction p s (T), which de- 
pends on temperature T, with p s (T = 0) = 1. Since we 
consider a 2D setup, a Berezinskii-Kosterlitz-Thouless 
transition towards a normal fluid is expected to occur at 
a finite temperature T KT = np s h 2 n/2m. In the opposite 



limit of strong interactions r,j > 1 the polar molecules 
are in a crystalline phase for temperatures T < T m with 
T m ps 0.09£>/a 3 ~ 0.018r d £' RiC , while for larger tempera- 
ture the crystal melts into a normal fluid via a first-order 
(classical) phase transition. The critical value T m for this 
melting transition has been obtained via molecular dy- 
namics simulations in the context of interfacial colloidal 
crystals in Ref^. Here, -Er ;C = ir 2 h 2 /2ma 2 is the crystal 
recoil energy, typically a few to tens of kHz. The config- 
uration with minimal energy is thus a triangular lattice 
with spacing a L = (4/3) ^a. Excitations of the crystal 
are acoustic phonons with Hamiltonian 

H c = '^huj q a' { q a q , (26) 
i 

where a q destroys a phonon of quasimomcntum q in the 
mode A. The characteristic Debye frequency is fiwo ~ 
l-6y/rdEji tC . At T = the static structure factor S di- 
verges at a reciprocal lattice vector K, and thus S(K)/N 
acts as an order parameter for the crystalline phase. 

In Refi^ we investigated the intermediate strongly in- 
teracting regime with > 1, and we determined the 
critical interaction strength r QM for the quantum phase 
transition between the superfluid and the crystal. In our 
analysis we used a recently developed PIMC-code based 
on the Worm algorithm^, which is an exact Monte-Carlo 
method for the determination of thermodynamic quan- 
tities in continuous space at small finite temperature. 
In Fig. [TOTd-e). the order parameters p s and S(K)/N 
are shown at a small temperature T = 0.014D/a 3 for 
different interaction strengths and particle numbers 
N = 36, 90. We find that p s exhibits a sudden drop 
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to zero for ra ~ 15, while at the same position S(K) 
strongly increases. In addition, during the Monte-Carlo 
simulations we observed that in a few occasions p s sud- 
denly jumped from to 1, and then returned to 0, in 
the interval f=a 15 — 20, which suggests a competition 
between the superfluid and crystalline phases. These re- 
sults indicate a superfluid to crystal phase transition at 



' QM 



18 ±4. 



(27) 



The step-like behavior of p s and S(K)/N is consistent 
with a hrst order phase transition, a result which 
has been confirmed in Refs . 51 ' 52 , We notice that the 
superfluid with r c i ~ 1 is strongly interacting, and 
in particular the density-density correlation function 
is quenched on lengths R < a, see Fig. [TOt ch This 
observation is consistent with the validity condition of 
the effective 2D interaction potential Eq. that two 

particles never approach each other at distances smaller 
than 

Having determined the low-tcmpcraturc phase- 
diagram, the remaining question is whether these 
phases, and in particular the crystalline phase emerging 
at strong dipolc-dipole interacions, are in fact accessible 
with polar molecules. This question is addressed in 
Fig. [TOf a), which is a contour plot of the interaction 
strength r d as a function of the induced dipole moment 
d = \J~D (in units of Debye) and of the mean interpar- 
ticle distance a (in /xm). The dimcnsionless quantity 
\J m/200amu depends on the mass m of the molecules 
(in atomic units), and it is of order one for characteristic 
molecules like SrO or RbCs. In the figure, stable 2D 
configurations for the molecules exist in the parameter 
region where the transverse (optical) trapping frequency 
oj = 2tt x 150Hz exceeds the dipolc-dipole interaction 
(hw > D/a 3 ), such that l ± = (12D/muj 2 ± ) 1/5 < a 
consistent with the stability discussion of Sect. IIII B 11 
[notice that l± ~ (D/huj±y^ 3 for realistic parameters]. 
The figure shows that for a given induced dipole d the 
ground-state of an ensemble of polar molecules is a 
crystal for mean interparticle distances l± < a < a max , 
where a max = <i 2 m/ft 2 7'QM corresponds to the distance 
at which the crystal melts into a superfluid. For SrO 
(RbCS) molecules with permanent dipole moment 
d = 8.9D (d = 1.25D), a min - 200nm(100nm), while 
Omax can be several ixm. Since for large enough inter- 
actions the melting temperature T m can be of order of 
several pK, the self-assembled crystalline phase should 
be accessible for reasonable experimental parameters 
using cold polar molecules. 

For what concerns the observability of the zero- 
temperature phases, Bragg scattering with optical 
light allows for probing the crystalline phase, while 
the detection of vortices can be used as a definitive 
signature of superfluidity. We notice that the 2D (quasi) 
condensate involves a fraction of the total density only, 
and therefore wc expect only small coherence peaks in a 
time of flight experiment. 
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FIG. 11: A dipolar crystal of polar molecules provides a peri- 
odic lattice V cp for extra atoms or molecules giving rise to a 
lattice model with hopping J and long-range interactions Vij 
(see text and Fig. [SJ. (a) A ID dipolar crystal with lattice 
spacing a provides a periodic potential for a second molecular 
species moving in a parallel tube at distance b (Configura- 
tion 1). (b) ID setup with atoms scattering from the dipolar 
lattice (Configuration 2). 



Finally, we notice that by adding an additional in- 
plane optical confinement, it is possible to realize strongly 
interacting ID phases which are analogous to the 2D 
crystals discussed abov o 53 ' 54 i 55 . For large enough in- 
teractions r 1, the phonon frequencies have the 
simple form hw q = {2/u 2 ) [12r d f q ] 1/2 £ RiC , with f q = 
X^>o ^ sin(qaj /2) 2 / j 5 . The Debye frequency is Rwd 



7r j a 



lA^/r2E^ c , while the classical melting tem- 



perature can be estimated to be of the order of T„ 
0.2r d E Kc /k B , see Ref^. 



B. Floating lattices of dipoles 

An interesting possibility offered by the realization of 
the self-assembled crystals discussed above is to utilize 
them as floating mesoscopic lattice potentials to trap 
extra-particles, which can be atoms or polar molecules 
of a different species. We show below that within an ex- 
perimentally accessible parameter regime extended Hub- 
bard models with tunable long-range phonon-mcdiatcd 
interactions describe the effective dynamics of the extra- 
particles dressed by the lattice phonons. 

The setups that we have in mind are shown in Figs. [5] 
and [Tl7a-b) where extra particles confined to a 2D crys- 
tal plane or a ID tube scatter from the periodic lattice 
potential • V cp (Rj — r). Here, r and Rj = R° + Uj 
are the coordinates of the particle and crystal molecule 
j, respectively, with R° the equilibrium positions and Uj 
small displacements. For particles being molecules, this 
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potential is given by the repulsive dipolc-dipole interac- 
tion y cp (Rj — r) = dpd c / | Rj — r | 3 with d p <C d c the 
induced dipole moment, and for atoms we assume that 
the interaction is modeled by a short range pseudopoten- 
tial proportional to an elastic scattering length a cp . In 
addition, extra molecules and atoms will interact accord- 
ing to dipolar, or short range interactions, respectively. 

We are interested in a situation where the extra par- 
ticles in the lattice are described by a single band Hub- 
bard Hamiltonian coupled to the acoustic phonons of the 
latticed 



H = j E c !c, + iE^ c H^ 



<i,3> 



(28) 



The first and second terms describe the nearest neighbor 
hopping of the extra particles with hopping amplitudes </, 
and interactions V, computed for each microscopic model 
by band-structure calculations for Uj = 0, respectively. 
The operators Cj (cj) are destruction (creation) operators 
of the particles. The third term is the phonon coupling 
obtained in lowest order in the displacement 

U 3 =iY J W^i c Nu q ) 1 ' 2 ^ q {a q + al g )e lq - R °, 



with 



M q = Vqq • ^(h^Nm^) 1 ^^ 



Here, £ 9 and N are the phonon polarization and the 
number of lattice molecules, respectively, while Vq is 
the Fourier transform of the particle-crystal interaction 
V cp , and /3 q = J c?r|wo(r)| 2 e lqr , with wq(t) the Wannier 
function of the lowest Bloch bancU^ 6 -. The validity of 
the single band Hubbard model requires J, V < A, and 
temperatures fcgT < A with A the separation to the 
first excited Bloch band. 



The Hubbard parameters of Eq. (|28p are of the or- 
der of magnitude of the recoil energy, J, V ~ Sr, C! and 
thus (much) smaller than the Debye frequency hwr> ~ 
ER,c^/rd, for rd 3> ti21. This separation of time scales 
J, V<^ huj£i, combined with the fact that the coupling to 
phonons is dominated by high frequencies ftio > J,V (see 
the discussion of M q below) is reminiscent of polarons 
as particles dressed by (optical) phonons, where the dy- 
namics is given by coherent and incoherent hopping on 
a lattic e) 118 ' 126 . This physical picture is brought out in a 
master equation treatment within a strong coupling per- 
turbation theory. The starting point is a Lang-Firsov 
transformation of the Hamiltonian H — > SHS^ with a 
density-dependent displacement 



S = exp 



E 

<!■.) 



Mi 



This eliminates the phonon coupling in the second line 
of Eq. (|28p in favor of a transformed kinetic energy term 



j e AcjXlXj 

<i,3> 

where the displacement operators 



Xj = exp 



M, 



can be interpreted as a lattice recoil of the dressed par- 
ticles in a hopping process. In addition, the bare inter- 
actions are renormalized according to 



V« = Vu 



r. 



(i) 



with V}p = -2 J2 q cos(q(R° - R°))M 2 /?ia; g , that is, the 
phonon couplings induce and modify off-site interactions. 
The onsite interaction is given by Vjj 



Vj.j - 2E p with 



E -V M ' 



the polaron self-energy or polaron shift. For J = the 
new Hamiltonian is diagonal and describes interacting 
polarons and independent phonons. The latter are vi- 
brations of the lattice molecules around new equilibrium 
positions with unchanged frequencies. A stable crystal 
requires the variance of the displacements Au around 
these new equilibrium positions to be small compared to 
a. 

A Born-Markov approximation with the phonons a fi- 
nite temperature heatbath with J, V <C hujo (see above), 
and the transformed kinetic energy 

j E - 

<i,j> 

as the system-bath interaction with {{XjXj)) the equi- 
librium bath average, provides the master equation for 
the reduced density operator of the dressed particles pt 
in Lindblad formi 2 ^ 



Pi 



[pt,H}+ E -fj~ ([bjS,Ptbis>] + [bis>,ptbjs]): 



3,1,8,8' 



(29) 



with bjs = Cj_^gCj. The effective system Hamiltonian 

6 = - J E ^ + \ E ViAcha, (so) 



<hj> 



'■;/ 



is of the extended Hubbard type, valid for J, Vij,E p < A. 
For E p > A, Eq.(3) should be derived via a multi-band 
approach. Coherent hopping of the dressed particles is 
described by 

J = J((Xjx j }} = Je X p(-S T ), 
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b/a 1 0.2 0.4 b/a 0.8 



FIG. 12: Configuration 1 [Fig. Illf a)]: Hubbard parameters 
for dp/rfc = 0.1 and m = m p . (a) Lattice depth Vo in 
units of -Er, p vs. b/a for rd = 50 and 500. Thick continu- 
ous lines: tight-binding region 4J < A. (b) Reduction factor 
So (dashed dotted lines) and polaron shift E p /J (solid lines), 
for 4J < A. (c) Continuous lines: phonon-mediated interac- 
tions Vjj +1 . Horizontal (dashed) lines: Vjj+i. (d) Contour 
plot of Vjj+i/2J (solid lines) as a function of b/a and rd- A 
single-band Hubbard model is valid left of the dashed region 
(4J, Vij < A), and right of the black region (E p < A). 

where 

St = E (j^) t 1 - cos(qa)](2n 9 (T) + 1) 

characterizes the strength of the particle-phonon interac- 
tions, and n q [T) is the thermal occupation at tempera- 
ture Ti26. 

The dissipative term in Lindblad form in Eq. (|29|) cor- 
responds to thermally activated incoherent hopping with 
rates T^'f , which can be made negligible for the energies 
of interest k B T < min(A, E p , k B T c ), see Refs . 118 i 129 . 
Corrections to Eq. ([501) proportional to J 2 are small 
relative to H provided J <C Ep^l& (also J <C hwo in 
l D 129 i 130 ). Thus, in the parameter regime of interest the 
dynamics of the dressed particles is described by the ex- 
tended Hubbard Hamiltonian H. In the following, we 
verify the existence of this parameter regime and we cal- 
culate the effective Hubbard parameters from the mi- 
croscopic model for the ID configuration of Fig. [TTT a). 
where extra-particles are polar molecules of a different 
species. An analogous calculation for the configuration 
of Fig. [TTT b) is reported in Refi^. 



In the configuration of Fig. [TIT a) molecules of a second 
species are trapped in a tube at a distance b from the 



crystal tube under ID trapping conditions. For crystal 
molecules fixed at the equilibrium positions with lattice 
spacing a, the extra particles feel a periodic potential 

V cp (x) = d c d p [b 2 + {x- ja) 2 ] , 
o 

which determines the bandstructure. The lattice depth 

d m P - 3b / a 

is shown in Fig. [T2T a) as a function of b/a, where the 
thick solid lines indicate the parameter regime 4 J < A, 
and £7r iP = ti 2 ir 2 /2m p a 2 . The potential is comb-like 
for b/a < 1/4, since the particles resolve the individ- 
ual molecules forming the crystal, while it is sinusoidal 
for b/a > 1/4. The strong dipole-dipole repulsion be- 
tween the extra particles acts as an effective hard-core 
constraint^ 9 -. We find that for 4J < A and d p <C d c the 
bare off-site interactions satisfy 

^•~4/(a|i-j|) 3 <A, 

which justifies a single-band approximation for the dy- 
namics of the extra-particles in the static potential. 
The particle-phonon coupling is 

with K-i the modified Bessel function of the second kind, 
and M q ~ Jq for q — > 0. In the regime of interest b/a < 1 
where the single-band approximation is valid (4 J, Vij < 
A), we find that M q is peaked at large q ~ 7r/a, so that 
the main contribution to the integrals in the definition 
of St and E p is indeed dominated by large frequencies 
hojq > J. Together with the separation of time-scales 
J, Vij *C fcfl, this is consistent with the picture of the 
system's dynamics as given by particles dressed by fast 
(optical) phonons, as discussed above. We notice that 
this so-called anti-adiabatic regime is generally hard to 
achieve in cold atomic setups^i. 

A plot of So as a function of b/a is shown in Fig.bLWb). 
We find the scaling 

So oc y/r2(d p /d c ) 2 , 

and within the regime of validity of the single band ap- 
proximation, So can be tuned from So «C 1 (J ~ J) to 
So 3> 1 (J J) corresponding to the large and small 
polaron limit, respectively The polaron shift E p gener- 
ally exceeds the bare hopping rate J, and in particular, 
E p > J for So > 1, see Fig. [IWb). Together with the 
condition Hlod 3> J this ensures that the corrections to 
Eq. (f50)) which are proportional to J 2 are indeed small, 
and thus Eq. (f5Uj) fully accounts for the coherent dynam- 
ics of the dressed particles. 

The extended Hubbard model corresponding to the 
configuration of FigfTlTa) is characterized by tunable off- 
site interactions, which arc a combination of the direct 
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dipolc-dipolc interactions between the extra-particles 
and of the phonon-mediated interactions v}^ . For b/a < 



decay slowly with 
2 and arc thus 



1 /4 we find that the interactions V i 
the inter-particle distance as ~ l/\i — j 
long-ranged. The sign of V i j is a function of the ratio 
b/a. Thus, depending on b/a the phonon-mediated in- 
teractions can enhance or reduce the direct dipolc-dipole 
repulsion of the extra particles. As an example, Fig. [T27 c) 

shows that the sign of the term Vjj +1 alternates between 
attractive and repulsive as a function of b/a, and that for 
small enough b/a the phonon-mediated interactions can 
become larger than the direct dipolc-dipole interactions. 

The effective Hubbard parameters Vjj+i and J are 
summarized in Fig. 112( d). which is a contour plot of 
Vj.j + i/2J as a function of and b/a. The ratio 
Vjj+i/2J increases by decreasing b/a or increasing Td, 
and can be much larger than one. This appearance of 
strong off-site interactions in the effective dynamics is 
a necessary ingredient for the realization of a variety 
of new quantum phase o 82 i 83 i 84 i 85 i 86 . As an example of 
the possible quantum phases that can be realized in this 
setup, at half filling, and considering nearest-neighbor 
interactions only, the particles in the configuration de- 
scribed above undergo a transition from a (Luttinger) liq- 
uid (Vis + i < 2 J) to a chargc-density-wave (Vi,i+i > 2J) 
as a function of b/a and r<j. Figure [T27d) shows that the 
parameter regime Vi.i+i ~ 2J, see RefJ^L, where this 
transition occurs can be satisfied for various choices of Td 
and b/a, e.g. for = 100 and b/a w 0.5. 



C. Three-body interactions 

As discussed in Sect. [Til it is of interest to design sys- 
tems where effective many-body interactions dominate 
over the two-body interactions, and determine the prop- 
erties of the groundstate. We here describe how an ef- 
fective low energy interaction potential V^ff of the form 
Eq. @ can be derived in the Born-Oppcnheimer approx- 
imation for £ polar molecules interacting via dipolc- 
dipole interactions by dressing low lying rotational states 
of each molecule with external static and microwave 
fields, in analogy to the discussion of Sect. IIII Bl 

We here focus on a setup with a static electric field 
E = Ee z along the z-axis, see Fig. [131 where the two 



states \g)i = |^ ,o)i and |e + ). t = \<j>i t +i)i with energies 2? fl 
and E e< ± are coupled by a circularly polarized microwave 
field propagating along the z-axis. The microwave tran- 
sition is characterized by the (blue) detuning A > and 
the Rabi frequency Sl/h. While the following discus- 
sion can be readily generalized to include the degenerate 
case^, here we assume that the degeneracy of the states 
|e_) = l^i.-i),: and |e+) is lifted, e.g., by an additional 
microwave field coupling the state |e_) near-resonantly 
to the next state manifold, see Fig. [T31 Then, the in- 
ternal structure of a single polar molecule reduces to a 



two- level system and is described as a spin- 1/2 particle 
via the identification of the state \g)i (|e+)$) as eigen- 
state of the spin operator Sf with positive (negative) 
eigenvalue. In the rotating frame and applying the ro- 
tating wave approximation, the Hamiltonian describing 
the internal dynamics of the polar molecule reduces to 



H, 



if a n 

2 I O -A 



(31) 



with the effective magnetic field h = (f2, 0, A) and the 
spin operator Sj = (Sf , Sf , Sf). The eigenstates of this 
Hamiltonian are denoted as = a\g)i + /3|e + ), and 
■ a|e+)j with energies iv^A 2 



n 2 /2. 



6B- 



2B- 



0- 




FIG. 13: Spectrum of a polar molecule, (a) Level structure 
for Ed/ B = 3: the circular polarized microwave field couples 
the ground state \g) with the excited state \e+) with Rabi 
frequency Q/h and detuning A. The excited state |e+) is 
characterized by a finite angular momentum J z \e+) = |e+). 
Applying a second microwave field with opposite polarization 
(red arrow) allows us to lift the degeneracy in the first ex- 
cited manifold by resonantly couple the state |e_) to the next 
manifold, (b) Internal excitation energies for a single polar 
molecule in a static electric field E = Ee z . 

For distances |r#| > (D/B) 1 / 3 with D = |(. 9 |d l |e + )| 2 
and d{ the dipole operator, the dipole-dipolc interaction 
Eq. ^ between two polar molecules can be mapped 
onto the effective spin interaction Hamiltonian H d = 
if™* + iPj hift . The first term describes an effective spin- 
spin interaction 



H? = -±Y l Dv(T ij )[SfSf 



SfS] 



')'- 



SfS]] , (32) 



>r.i 



where r]± = r\ g ± r\ e is determined by the induced dipole 
moments rj g = dEE g /\fT) and r\ e = dEE e ^+/VD. The 
anisotropic behavior of the dipole-dipole interaction is 
accounted for by v(r) = (1 — 3 cos 2 -d)/r 3 with ■& the angle 
between r and the z-axis. In addition, the asymmetry 
of the induced dipole moments gives rise to a position 
dependent renormalization of the effective magnetic field 
and a energy shift 



2 fe 



Dv(vij 



V-V+ 



n 2 

Sf + ^± 



(33) 
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Within the Born-Oppcnhcimcr approximation, an an- 
alytic expression for the effective interaction V^r^r,}) 
between two polar molecules each prepared in the state 
can be derived in second-order perturbation theory 
in the dipole-dipole interaction Vdd(r)/h as 



K ff ({rJ)-£ (1) ({rJ) +E^({r l }) 



(34) 



where £>/(et 3 |h|) = (i?o/a) 3 <C 1 is the (small) parameter 
controlling the perturbative expansion, a is the charac- 
teristic length scale of the interparticle separation and 
i?o = (-D/a/A 2 + fi 2 ) 1 / 3 is a Condon point, analogous to 
that discussed in Sect. lIIIB2l The energy shift 



(35) 

gives rise to a dipole-dipole interaction between the par- 
ticles, while the term 



E {2) ({r<}) 



E 



\M\ 



VA 2 + fi 2 



D 2 v{Y ik )v(v ]k ) 



E 

l<j 



\N\' 



2 V / A 2 + n 2 



[Dv (r, 



(36) 



corresponds to a correction to the two-particle interaction 
potential and an additional three-body interaction. The 
matrix elements M and N take the form 

M = ap [(a 2 Vg + p 2 Ve ) (rj e - t} g ) - (« 2 - /3 2 )/2] , 

N = a 2 f3 2 [(r/ e - Vg f + 1 

Therefore, the effective interaction potential V c{t up to 
second order in (Ro/a) 3 reduces to the form in Eq. ([2|) 
with the two-particle interaction potential 



V(r)=\ 1 Du(r) + X 2 DRl [v (r)] 5 
and the three-body interaction 
W(ri,r 2 ,r 3 ) = j 2 RqD [v(r l2 )i/(r 13 ) 



(37) 



(38) 

+ K r 12)K r 23) + K r 13)K r 23)] ■ 

The dimensionless coupling parameters are Ai = 
(a 2 Vg + /3 2 Ve ) 2 - a 2 /3 2 7 A 2 = 2|M| 2 + \N\ 2 /2, and 72 = 
2|M| 2 . These parameters can be tuned via the strength 
of the electric field Ed/B and the ratio between the Rabi 
frequency and the detuning, fi/A, see Fig. [14] Of spe- 
cial interest are the values of the external fields where 
the leading two-particle interaction vanishes, i.e. Ai = 0. 
Then, the interaction is dominated by the second order 
contribution with A 2 and 72, which includes the three- 
body interaction, see Fig. I14f d). while a small deviation 
away from the line Ai = allows us to change the char- 
acter of the two-particle interaction. Note, that a n-body 
interaction term (n > 4) appears in (n — l)-th order per- 
turbation theory in the small parameter (i?o/a) 3 . There- 
fore, the contribution of these terms is suppressed and 
can be safely ignored. 
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FIG. 14: Parameters of the effective interaction potential. 
(a)-(c): Strength of the interaction parameters Ai, A2, and 72 
as a function of the external fields Ed/B and fi/ A. The lead- 
ing dipole-dipole interaction vanishes for Ai = [dashed line 
in (b) and (c)], and the second order contributions dominate 
the interaction, (d) Strength of A2 (dashed line) and 72 (solid 
line) along the line in parameter space with Ai = 0. 



The perturbative expansion requires that a 3> -Ro- 
We notice that particles can be confined to interparticle 
distances larger than Rq by combining repulsive dipole- 
dipole interactions with a strong (optical) transverse con- 
finement ui±, in analogy to Sect. IIII B II Then, for a 



repulsive two-particle potential with Ai 



> 



-\ 2 {Ro/aY 



and for Hlj±_ > D/Rq, two-particles reach distances 
|rj — r j I < Ro at an exponentially small rate T ~ 
(h/ma 2 )exp(-2S B /h), with S E /H ~ yjDm/Roh 2 . This 
exponential suppression ensures the stability of the colli- 
sional setup for the duration of an experiment. 

The low-energy many-body theory now follows by com- 
bining the kinetic energy of the polar molecules with the 
effective interaction V o{t within the Born-Oppcnheimcr 
approximation and the external trapping potentials V T 



^=E^+^«(M)+E^H- 



(39) 



This Hamiltonian is independent of the statistics of 
the particles and therefore it is valid for bosonic and 
fermionic polar molecules. 

Applying an optical lattice provides a periodic struc- 
ture for the polar molecules described by the Hamiltonian 
Eq. (|3"5|) . In the limit of a deep lattice, a standard ex- 
pansion of the field operators i/j^(r) = J2i w ( r ~ 
in the second-quantized expression of Eq. (|39[) in terms 
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of lowest-band Wannier functions w(r) and particle cre- 
ation operators &ti2L leads to the realization of the Hub- 
bard model of Eq. ©, characterized by strong nearest- 
neighbor interactions^. We notice that the particles are 
treated as hard-core because of the constraint o» JJo, 
The interaction parameters Uij and Vyfc in Eq. ((9]) derive 
from the effective interaction V({ri}), and in the limit of 
well-localized Wannier functions reduce to 



U« = U - 



R, 



and 



w, 



Wo 



IRi-RjPIRi-Rit 



+ perm 



(40) 



(41) 



respectively, with [/o = XiD/a 3 , U\ = A2-D-R0/ 'a 6 , and 
Wo = I2DR0I 'a 6 . The dominant contributions and 
strengths of the three-body terms in different lattice ge- 
ometries are shown in Fig. H£b). For LiCs with a perma- 
nent dipole moment d = 6.3Debye trapped in an optical 
lattice with spacing a « 500nm, the leading dipolc-dipole 
interaction can give rise to very strong nearest-neighbor 
interactions with Uq ~ 55i? kin , and E kiri = h 2 /ma 2 . On 
the other hand, tuning the parameters via the exter- 
nal fields to Ai = the characteristic energy scale for 
the three-body interaction becomes Wo « (Ro/a) 3 E kia . 
Then, controlling the hopping energy J via the strength 
of the optical lattice allows to enter the regime with dom- 
inant three-body interactions. For particles been bosons, 
an analytic calculation has suggested that the ground- 
state phase diagram of Eq. © with Uij = in ID is char- 
acterized by the presence of valence bond states at spe- 
cific rational fillings of the lattice, charge-density waves 
and supcrfluid phases^. 



D. Lattice Spin models 

Cold gases of polar molecules allow to construct in a 
natural way a complete toolbox for any permutation sym- 
metric two spin-1/2 (qubit) interaction, using techniques 
of interaction engineering similar to those discussed in 
the previous sections. 

The system we have in mind is comprised of het- 
eronuclear molecules with 2 T, 1 / 2 ground electronic states, 
corresponding for example to alkalinc-carth monohalo- 
genides with a single electron outside a closed shell. We 
adopt a model molecule where the rotational excitations 
are described by the Hamiltonian 



H m = BN 2 + 7 N • S, 



(42) 



with N the dimensionless orbital angular momentum of 
the nuclei, and S the dimensionless electronic spin (as- 
sumed to be S — 1/2 in the following). Here B denotes 
the rotational constant and 7 is the spin-rotation cou- 
pling constant, where a typical B is a few tens of GHz, 



and 7 in the hundred MHz regime. The coupled basis of 
a single molecule i corresponding to the eigenbasis of H* n 
is {\Ni,Si, JjjMjJ} where J,; = N, + Sj with eigenval- 
ues E(N = 0,1/2,1/2) = 0, £(1,1/2, 1/2) =2B- 7 , and 
£(1,1/2,3/2) = 2B + 7/2. 

The Hamiltonian describing the internal and external 
dynamics of a pair of molecules trapped in wells of an 
optical lattice is denoted by H = H m + H cx . The in- 
teraction describing the internal degrees of freedom is 
H m = i?dd + Si=i -^m; where H^d is the dipole-dipole 
interaction. The Hamiltonian describing the external, or 
motional, degrees of freedom is H cx = Y^i=i P?/(2n<) + 
Vi(xi— Xj), where Pi is the momentum of molecule i with 
mass m, and the potential generated by the optical lattice 
Vi(x — xj) describes an external confinement of molecule 
i about a local minimum Xj with ID rms width zq. We 
assume isotropic traps that are approximately harmonic 
near the trap minimum with a vibrational spacing hio osc . 
Furthermore, we assume that the molecules can be pre- 
pared in the motional ground state of each local potential 
using dissipative electromagnetic pumpingi22,. It is con- 
venient to define the quantization axis z along the axis 
connecting the two molecules, X2 — xi = Azz with Az 
corresponding to a multiple of the lattice spacing. 

The ground subspace of each molecule is isomor- 
phic to a spin 1/2 particle. Our goal is to obtain 
an effective spin-spin interaction between two neigh- 
boring molecules. Static spin-spin interactions due 
to spin-rotation and dipole-dipole couplings do ex- 
ist but are very small in our model: -ff v dw( r ) = 
-(d 4 /2Br 6 ) [l + ( 7 /4B) 2 (1 + 4Si • S 2 /3 - 2SfS$) . 

The first term is the familiar van-der-Waals 1/r 6 in- 
teraction, while the spin dependent piece is strongly 
suppressed as 7/4B sa 10~ 3 <C 1. However, dipole-dipole 
coupled excited states can be dynamically mixed using 
a microwave field. 

The molecules are assumed to be trapped with a sep- 
aration Az ~ r 7 = (2d 2 /7) 1 / 3 , where the dipole dipole 
interaction is d 2 /r 3 — 7/2. In this regime the rotation 
of the molecules is strongly coupled to the spin and the 
excited states are described by Hunds case (c) states in 
analogy to the dipole-dipole coupled excited electronic 
states of two atoms with fine-structure. The ground 
states are essentially spin independent. In the subspace 
of one rotational quantum (N\ + N2 = 1), there are 24 
eigenstates of H ln which are linear superpositions of two 
electron spin states and properly symmetrized rotational 
states of the two molecules. There are several symme- 
tries that reduce _ff; n to block diagonal form. First, 
i/dd, conserves the quantum number Y = Mn + Ms 
where M N = M Nl + M N2 and M s = M Sl + M S2 
are the total rotational and spin projections along the 
intermolecular axis. Second, parity, defined as the 
interchange of the two molecules followed by parity 
though the center of each molecule, is conserved. The 
a = ±1 eigenvalues of parity arc conventionally denoted 
g(u) for gcradc(ungcradc). Finally, there is a symme- 
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try associated with reflection R of all electronic and 
rotational coordinates through a plane containing the 
intermolecular axis. For \Y\ > all cigenstates are even 
under R but for states with zero angular momentum 
projection there are ±1 eigenstates of R. The 16 distinct 
eigenvalues correspond to degenerate subspaces labeled 
|y|^(J) with J indicating the quantum number in the 
r — > oo asymptotic manifold (N = 0, J = 1/2; N = 1, J). 
Remarkably, the eigenvalues and eigenstates can be 
computed analytically yielding the Movrc-Pichler poten- 
tials^ plotted in Fig.rjj^a). 

In order to induce strong dipole-dipole coupling we 
introduce a microwave field E(x,t)ep with a frequency 
lof and Rabi- frequency tuned near resonance with the 
N = -> N = 1 transition. 

The effective Hamiltonian acting on the lowest-energy 
states is obtained in second order perturbation theory as 



hf A(r) 



(g / |g mf |A(r))(A(r)|g mf | gi 
hoj F -E(X(r)) 



\9f)(9i 



N, 



(43) 
= 



where {\gi), |<7/)} are ground states with N\ 
and {|A(r)}} are excited eigenstates of H ln with N\ + 
N2 = 1 and with excitation energies {E(A(r))}. The 
reduced interaction in the subspace of the spin degrees 
of freedom is then obtained by tracing over the mo- 
tional degrees of freedom. For molecules trapped in the 
ground motional states of isotropic harmonic wells with 
rms width zq the wave function is separable in center of 
mass and relative coordinates, and the effective spin-spin 
Hamiltonian is H spm = (H cS (r)) rcl . 

The Hamiltonian in Eq. ([43]) is guaranteed to yield 
some entangling interaction for appropriate choice of field 
parameters but it is desirable to have a systematic way to 
design a spin-spin interaction. The model presented here 
possesses sufficient structure to achieve this essentially 
analytically. The effective Hamiltonian on molecules 1 
and 2 induced by a microwave field is 



H eS (r) 



h\n\ 



°iA a Ar)° 

a, (3=0 



(44) 



where {u a }^ =0 = {1, cr x , <r v , a z } and A is a real symmet- 
ric tensor. 

Equation (|44[) describes a generic permutation sym- 
metric two qubit Hamiltonian. The components Aq iS 
describe a pseudo magnetic field which acts locally on 
each spin and the components A s> t describe two qubit 
coupling. The pseudo magnetic field is zero if the mi- 
crowave field is linearly polarized but a real magnetic 
field could be used to tune local interactions and, given 
a large enough gradient, could break the permutation in- 
variance of H spin . 



For a given field polarization, tuning the frequency 
near an excited state induces a particular spin pattern 
on the ground states. These patterns change as the fre- 
quency is tuned though multiple resonances at a fixed 



TABLE I: Some spin patterns that result from Eq. (|44|) . The 
field polarization is given with respect to the intermolecular 
axis z and the frequency ujf is chosen to be near resonant 
with the indicated excited state potential at the internuclear 
separation Az. The sign of the interaction will depend on 
whether the frequency is tuned above or below resonance. 
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intermolecular separation. In Table U it is shown how 
to simulate the Ising and Hcisenbcrg interactions in this 
way. Using several fields that are sufficiently separated 
in frequency, the resulting effective interactions arc addi- 
tive creating a spin texture on the ground states. The 
anisotropic spin model Hxyz = X x cr x o' x + \<J v a v + 
X z a z a z can be simulated using three fields: one polarized 
along i tuned to 0+(3/2), one polarized along y tuned to 
0~(3/2) and one polarized along y tuned to 0+ ( 1/2). The 
strengths A^ can be tuned by adjusting the Rabi frequen- 
cies and detunings of the three fields. Using an external 
magnetic field and six microwave fields with, for exam- 
ple, frequencies and polarizations corresponding to the 
last six spin patterns in Table U arbitrary permutation 
symmetric two qubit interaction are possible. 

The Kitaev model of Eq. (jllj) (Spin model II) can be 
obtained in the following way. Consider a system of four 
molecules connected by three length b edges forming an 
orthogonal triad in space. There are several different 
microwave field configurations that can be used to realize 
the interaction H^ n along the links. One choice is to 
use two microwave fields polarized along z, one tuned 
near resonance with a \ g potential and one near a l u 
potential. A realization of model II using a different set 
of 3 microwave fields is shown in Fig.fToTb'). The obtained 
interaction is close to ideal with small residual coupling 
to next nearest neighbors. 
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FIG. 15: (a) Movre-Pichler potentials for a pair of molecules as a function of their separation r: The potentials E(gi(r)) for 
the 4 ground-state (dashed lines) and the potentials E(X(r)) for the first 24 excited states (solid lines). The symmetries \Y\f 
of the corresponding excited manifolds are indicated, as are the asymptotic manifolds (Ni, Ji\ Nj , Jj). (b) Implementation of 
spin model H^ n . Shown is the spatial configuration of 12 polar molecules trapped by two parallel triangular lattices (indicated 
by shaded planes) with separation normal to the plane of Az/y/3 and in plane relative lattice shift of AzW2/3. Nearest 
neighbors are separated by b = Az and next nearest neighbor couplings are at y/2b. The graph vertices represent spins and 
the edges correspond to pairwise spin couplings. The edge color indicates the nature of the dominant pairwise coupling for 
that edge (blue= o z o z , red= o v o v , green= a x o x , black= "other" ) . For nearest neighbor couplings, the edge width indicates 
the relative strength of the absolute value of the coupling. For this implementation, the nearest neighbor separation is b — r 7 . 
Three fields all polarized along z were used to generate the effective spin-spin interaction with frequencies and intensities 
optimized to approximate the ideal model H^ n . The field detunings at the nearest neighbor spacing are: hui — _B(l 9 (l/2)) = 
-0.057/2, hu 2 - E(Qg (1/2)) = 0.057/2, hw 3 - E(2 g (3/2)) = O.lOy/2 and the amplitudes are = 4|fi 2 | = |n 3 | = 0.017/ft. 
For 7 = 40MHz this generates effective coupling strengths J z — — 100kHz and J± = — 0.4J Z . The magnitude of residual nearest 
neighbor couplings are less than 0.04|J Z | along x and y-links and less than 0.003|J Z | along z-links. The size of longer range 
couplings Ji r are indicated by edge line style (dashed: |Ji r | < 0.01 |J Z |, dotted: |Ji r | < 10~ 3 |J Z |). Treating pairs of spins on 
z-links as a single effective spin in the low energy sector, the model approximates Kitaev's 4-local Hamiltoniani 1 ^ on a square 
grid (shown here are one plaquette on the square lattice and a neighbor plaquette on the dual lattice) with an effective coupling 
strength J cB = -( J_l/J z ) 4 | J*|/16 « 167Hz. 
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